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ABSTRACT 


The Global Positioning System (GPS) is revolutionizing the science of navigation. Never before has 
there been a system that could provide real time, world wide, continuous coverage with the such precision. 
Yet, the accuracy achievable with GPS alone is not sufficient to achieve autonomous flight of an Unmanned 
Aerial Vehicle. However, when integrated with an Inertial Navigation System and other non-inertial sensors 
using a Kalman Filter, GPS supplies the critical positioning information to permit such an achievement. 

This thesis presents the selection and integration of a GPS receiver using Differential GPS (DGPS) 
in support of a UAV autonomous flight project. Contemporary electronic navigation systems are surveyed, 
GPS operation is reviewed, and a Motorola PVT-6 GPS receiver selected. Using the Motorola Proprietary 
Binary Format protocol, several software drivers Were written in C to interface the information to an Intel 
80486DX CPU using the RS-232 serial communication standard. Finally, an examination 1s made to 
determine the maximum reacquisition time, the DGPS accuracies achievable and the effects of pseudorange 


correction latency on DGPS accuracy. 
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I. INTRODUCTION 


The last of a 24 satellite constellation was launched into 
orbit in June 1993, and with it, "the science of navigation is 
poised to take a great leap forward" (Carey, 1993, p. Bl). The 
Global Positioning System (GPS) revolutionized electronic 
navigation by providing worldwide, real-time, continuous 
coverage with three dimensional accuracies at least an order 
of magnitude better than any other system in existence. Even 
with these capabilities, it is not the panacea some make it 
out to be. Absolute accuracies (1.e., accuracies using a 
Single receiver) better than ten meters within two standard 
deviations cannot be achieved. 

Differential GPS (DGPS) shows great promise by further 
improving accuracies to within one to three meters. DGPS uses 
two GPS receivers where one 1S placed ina known location and 
the other is allowed to move. By computing the difference 
between the stationary receiver’s actual and GPS determined 
locations, and then transmitting these differences to the 
other GPS receiver, the other GPS receiver can correct its 
output to obtain the improved position information. This 
improvement would be sufficient for an aircraft to achieve an 
autonomous landing with GPS alone were it not for the 
unreliability of the GPS signals in high speed, critical 


maneuvering evolutions. It is not uncommon for a GPS receiver 


to lose track of one or more of its satellites for a several 
seconds, and thereby, be unable to provide positioning 
information. While not a catastrophe for a slow moving surface 
vehicle, an outage of even a few seconds could prove fatal for 
an aircraft. 

In an attempt to achieve the accuracy and reliability 
required for autonomous flight in a navigation package light 
enough for a small unmanned aerial vehicle (UAV), this 
research effort is part of a project where the GPS information 
is combined with the outputs of an Inertial Measurement Unit 
(IMU) and with data from other non-inertial navigation sensors 
such as airspeed and altitude indicators in a Kalman Filter. 
A contribution is made by selecting the best GPS receiver, 
defining its method of employment and writing the necessary 
software drivers to smoothly interface them with a Central 
Processing Unit (CPU). 

In the following’ chapter, contemporary electronic 
navigation systems are surveyed, discussing their advantages 
and disadvantages. Realizing the obvious benefits of GPS, 
Chapter III examines the User Segment of GPS emphasizing those 
details necessary to intelligently select and use the GPS 
receiver in a manner that best supports the various project 
requirements. Chapter IV brings it all together by detailing 
the receiver selection and capabilities as well as the 
criteria for its use. It also describes the hardware and 


software necessary to interface the receivers with the CPU. 


Three tests were conducted to confirm expected operational 
Capabilities. They are described and the results from them are 
presented in Chapter V. The appendices include receiver 
Specifications, the ephemeris algorithm and the software 
drivers written in C code as well as other information to aid 


the user and the programmer. 


IIT. SURVEY OF ELECTRONIC NAVIGATION METHODS 


The invention of wireless communications in the early 
twentieth century ushered in a new era in navigation allowing 
seafarers to safely pilot their ships in the absence of sight 
and sound for the first time in the history of mankind. The 
first use of radio waves required a vessel to transmit a 
continuous wave at some frequency. On shore, receivers with a 
rudimentary direction finding capability triangulated the 
ship’s position and transmitted it back to the ship. The 
position was not very accurate by today’s standards but served 
its purpose well in regions of the world where the sky was 
commonly obscured by clouds or fog when far from land. 

The next significant advance came with the invention of 
highly accurate chronometers. First, the quartz clock, and 
then, the atomic time standards with mind boggling stabilities 
providing accuracies better than one second in 30,000 years or 
one part in a trillion (107) (Logsdon, 1992, p.156). Such 
precision is required to allow the precise measurement of 
frequency. in adadaeaon- Transit (SatNav) and Global 
Positioning System (GPS) partially determine position by 
measuring the amount of time for the signal to travel from the 
satellite to the receiver. Radio waves travelling at the speed 


of light will advance one foot in one billionth of a second 


mavcon, 1990, p.336)2 .-It is not. difficult to see the 
necessity for such accuracy in time Keeping. 

In order to appreciate the quantum improvement in the 
science of navigation that GPS has achieved, it is worthwhile 
to review the major electronic navigation methods in use 
today. A brief review of some principles and characteristics 
of electromagnetic waves is given, followed by an individual 
description of each of the major navigation systems. While an 
effort has been made to succinctly describe the theory upon 
which the system operates, the important point is to notice 


the range, coverage and accuracy for each method. 


A.- PRINCIPLES AND CHARACTERISTICS OF ELECTROMAGNETIC WAVES 
Radio waves in a vacuum theoretically travel at the speed 


of light. Wavelength is related to frequency by Equation (1): 
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where A is the wavelength in meters and f is the frequency in 
Hertz. Figure 1 demonstrates this relation for’ the 
electromagnetic spectrum. Information may be carried on a 
radio wave through one or more modifications to its shape. 
These modifications include changing the amplitude, frequency 
or phase of the radio wave and are appropriately called 
amplitude, frequency and phase modulation, respectively. In 


addition, a wave may be intermittently transmitted or pulsed. 


This is called pulsed modulation. Figure 2 shows’ three 


representative examples of modulation techniques. 
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Figure 1. Wavelength Related To Frequency In The 
Electromagnetic Spectrum. (Hobbs, 1981, p. 217) 


Radio waves, like light, experience such phenomena as 
reflection, absorption, refraction, diffraction and 
interference. With respect to the earth’s environment, the two 
most influential features on the distance and path radio waves 


travel are the earth’s surface and the ionosphere. The 
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Figure 2. Modulation Techniques Commonly Used. 
(Hebbs, 1981, p.219) 


lonosphere has four layers and exists at a height of 
approximately 60 to 400 Km above the earth’s surface. This 
highly charged and fluctuating segment of the atmosphere not 


only absorbs radio wave energy but also reflects it (Beck, 


1971, p.95). Waves that are reflected by the ionosphere are 
called "sky waves". The ionosphere is affected by diurnal and 
seasonal fluctuations and by the level of sunspot activity 
which in turn affects the performance characteristics of sky 
waves. Sky wave phenomena is also greatly influenced by the 
frequency, angle of incidence, and the height and density of 
the various layers of the ionosphere. The sky wave frequencies 
most commonly observed are in the 300 KHz to 30 MHz band 
(Medium and High Frequency bands). Because of the higher 
resistance of the earth’s crust compared to the atmosphere, 
radio waves transmitted parallel to the earth’s surface are 
slowed causing the wave to bend. The degree of bending is 
inversely proportional to the frequency. Thus lower 
frequencies will travel over the horizon in a direct path. 
These radio waves are called ground waves. As it might be 
expected, a receiver might travel along the surface of the 
earth from the transmitter origin and receive the ground wave, 
enter into an area shadowed by the earth’s curvature with no 
Signal, and then pick up the reflected "one-hop" sky wave even 
PUbeEher Cue. Figure 3 illustrates this. Furthermore, when 
conditions are right, one-hop sky waves may experience 
reflection by the earth’s surface and subsequent ionospheric 
reflection for a "two-hop" sky wave. Navigation systems such 
as Loran-C and Omega use reflected sky waves (when they exist) 


to extend their effective range. (Logsdon, 1992, pp. 95-98) 





Figure 3. Sky Wave And Ground Wave Patterns. Ray 1 Is 
Refracted Through The Ionosphere, Ray 2 Is A One-Hop 
Skywave, Ray 3 Is A Two-Hop Skywave And Ray 4 Is A Ground 
Wave. (Hobbs, 1981, p. 225) 


Finally, several major navigation systems exploit radio 
wave phase measurement in order to determine two dimensional 
position. The principle of operation uses two synchronized 
transmitters located a fixed distance apart. The line between 
them 1S called the baseline. As a result, at half-wavelength 
intervals the phase difference between the two signals is 
zero. In this type of system, the half-wavelength segments are 
normally called lanes. As the distance from the baseline 
increases the lines of equal phase or "isophase lines" radiate 
outward from the foci of the two transmitters in hyperbolic 
fashion as shown in Figure 4. Thus, methods employing this 
technique are called hyperbolic radio navigation systems. If 
the two transmitters were to transmit a pulse instead of a 
continuous wave, then by measuring the time of arrival of each 
synchronized pulse, one’s position might also be determined. 
The points where the times of arrival of transmitted pulses 


are the same describe a hyperbolic locus. This type of method 


1s also considered a hyperbolic radio navigation system. With 
the exception of radio beacons, all contemporary terrestrial 
radio navigation systems use some sort of hyperbolic 


navigation method. (Bowditch, 1984, pp. 1015-1016) 
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Figure 4. Construction Of A Hyperbolic Interference 
Navigation Pattern. (Hobbs, 1981, p. 227) 
B. RADIO BEACONS 

Radio beacons are continuous wave transmitters modulated 
with a unigue Morse code identification. They are usually low 
power transmissions in the low and mid frequency bands with a 


range of less than 200 miles. Most cannot be received beyond 


10 


20 miles (Hobbs, 1981, p. 231). They operate using shipboard 
radio direction finding (RDF) equipment where two or more 
beacons are used to obtain a fix. The accuracy is highly 
dependent on the sensitivity of the RDF equipment as well as 
the range from the beacons. As Figure 5 illustrates, the 
geometry of the radio beacons relative to the ship may also 
play a Significant role in the accuracy. Finally, the precise 
positions of the beacon transmitters must be known or the line 
of bearing is useless. In other words, the _ position 
information cannot be determined in terms of latitude or 


longitude until after it is plotted. 


C. LORAN-C 

Loran-C is the descendent of the first attempt at a 
hyperbolic navigation system developed during World War II for 
both ships and aircraft (Wilkes, 1987, pp. 31-34). From Figure 
6, 1t can be seen that it covers almost the entire northern 
hemisphere with only a very small area covered below the 
equator and then only by sky waves when they exist. 

Loran-C operates in the low frequency band at 100 KHz with 
a 10 KHz bandwidth. Within geographic areas, there is one 
master station and two or more secondary stations. All 
stations sequentially transmit a pulse modulated continuous 
wave in groups of eight pulses separated by one microsecond 
with the exception of the master station which transmits nine 


pulses. The ninth pulse is separated from the first eight by 


11 


. = . = € ¢ 
* er ssen she « “70m meee oe 


PLOTTING INTERSECTING RADIO SEARINGS 


Ship's most probable position... ~—@ 


Other probable positions . 





Figure 5. Intersection Of Radio Beacon Signals Showing Error 
Spread As Range From The Source Increases. (Hobbs, 1981, 
Din235} | 


two microseconds and is used primarily for visual 
identification when viewing the pulses on an oscilloscope. 


This can be seen in Figure 7. Each group of master and 
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Figure 6. Areas Of The World Covered By Loran-c. In 
Particular, Note That The Coverage Is Almost Completely In 
The Northern Hemisphere. (Hobbs, 1981, p. 237) 
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secondary stations is called a chain. Currently, there exist 
13 chains throughout the world. (USNO, 1972, p. 7-5) (Hobbs, 
1981, pp. 236-239) 

Each master and secondary station transmits according to 
a predetermined fixed order synchronized using precise atomic 
time standards. Within each receiver, the group of eight 
pulses are integrated to form one pulse of 320 microseconds. 


The time difference of arrival is then measured between the 
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Figure 7. Loran-C Pulse Sequence For A Four Station Chain. 
(Hobbs 7 193i p.e2c7) 


master and selected secondary station pulse from which may be 
plotted a line of position on Loran-C charts. As shown in 
Figure 8, Loran-C charts have the number coded hyperbolic 
lines overlaid on charts and to enable plotting. Chain and 
station identification are achieved by measuring the pulse 
repetition rate (PRR) and the pulse repetition interval (PRI). 
(USNO, 1972, p. 7-5) (Hobbs, 1981, pp. 236-239) 

Accuracy of the system varies from about 700 ft near the 
baseline to about 2000 ft near the extreme range of the 
system. Due to the low frequency, baseline distances range 
from 1000 to 1500 miles so that coverage within a chain 
includes a great deal of area. Similar to the use of radio 
beacons, Loran-C charts and tables are required for the signal 


information to be meaningful. 
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Figure 8. A Portion Of A Loran-C Chart Covered By The 
Northeast U.S. Chain. (GRI 9960) (Hobbs, 1981, 242) 


D. DECCA NAVIGATOR SYSTEM 

Decca was conceived by an American and developed by the 
British Admiralty Signals Establishment during World War II 
for use in mine sweeping and the Normandy invasion forces in 
1944 (Beck, 1971, p. 65). Today, it 1S operated by the Racal- 
Decca Navigator Company in only a few but very heavily 
traveled areas in the world as shown in Figure 9. It operates 
using the principles of hyperbolic radio navigation but on a 
much smaller scale than Loran-C. 

Each separate area of coverage has a master station and up 


to three slave stations arranged ina star pattern around the 
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Figure 9. Areas Of The World Covered By Decca. (Beck, 1971, 
p. 64) 


master station and identified by a color designator of green, 
red or purple. Each station within master/slave chain 
transmits an integer multiple of an unmodulated, continuous 
wave base frequency (f) in the range of 14.00 KHz to 14.33 
KHz. The master station transmits at 6f, the purple station at 
5f, the red station at 8f and the green station at 9f. Each 
master/slave pair is phase locked along the baseline. A Decca 
receiver consists of four receivers, one for each of the 
frequencies transmitted by the master and slave stations. 
(Beck, 1971, pp. 65-67) 

Within the receiver, each master/color frequency is 
integer multiplied to get another unigue multiple of the base 


frequency which is the first common multiple of either of the 
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two. For instance, the master station frequency of 6f when 
multiplied by four and the red station frequency of 8f when 
multiplied by three produce a new harmonic of 24f for the 
Master/red station pair. Each half wave-length of this new 
harmonic is considered a "Decca lane", and by comparing phases 
of the two signals, one can determine position within the lane 


as Figure 10 illustrates. 
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Figure 10. Phase Comparison Of Decca Master/Red-Slave 
Signals Along The Base Line. (Hobbs, 1981, p. 251) 


Lane ambiguities are resolved once each twenty seconds 
when all four stations transmit the fundamental frequency f 
whose half-wavelength precisely encompasses the common 
multiples (18f for the master/green pair, 24f for the 
master/red pair and 30f for the master/purple pair). By using 
phase comparison of the received signal in each of these 
"Decca zones" as shown in Figure 11, the exact lane can be 


determined. In a like manner, zone ambiguities are resolved by 
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transmitting on a frequency of 8.2f also every 20 seconds. 
Called an orange frequency, it forms a coarse hyperbolic 
pattern within which each 360° difference cycle encompasses 


five zones. (Hobbs, 1981, pp. 252-253) 
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Figure 11. Phase Comparison Of Decca Master/Green-Slave Lane 
Identification Signals. (Hobbs, 1981, p. 252) 


The effective range of each chain is approximately 240 
miles from the master station. At 100 miles, the accuracy is 
approximately 30 yds by day and 100 yds by night. Again, this 
system requires the Decca charts and receiver in order to 


provide positioning data. (Bowditch, 1984, pp.1058-1059) 


E. CONSOL 

Consol was developed in 1945 from a German radio beacon 
system called Sonne. It is a hyperbolic radio navigation 
system with a very short baseline such that a "degenerated" 


hyperbolic pattern is formed. The curved portions of the 
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hyperbolic pattern are not used, and so it acts as a long 
range radio beacon. It is used mostly in Europe and the 
Mediterranean. Another system very Similar to Consol, called 
Cosolan, iS used in the United States. Both have the advantage 
that no special equipment iS required other than a low 
frequency receiver. The operator only needs to listen to the 
number of dots or dashes to know where within a particular 
sector the receiver is located. (Beck, 1971, 113-116) 

Both systems use low frequency radio waves in the range 
from 190 KHz to 350 KHz. Each station consists of three radio 
towers spaced three wavelengths apart. One of the towers 
transmits a continuous wave frequency while the other two 
transmit waves that undergo 180° phase shifts with respect to 
the CW wave during a keying cycle. All the signals are 
modulated with dots and dashes. This arrangement produces a 
series of sectors 10° to 15° wide as shown in Figure 12. The 
phase shift of the combined signal results in a variable 
number of dots and dashes being heard depending on the radial 
position within the sector. (Beck, 1971, 113-116) 

Consol range can be up to 1500 miles but not less than 25 
miles due to the compressed width between sectors close to the 
transmitters. Keeping in mind that only radial lines of 
position may be obtained from a Consol station, accuracies of 
less than a degree may be obtained at the maximum range 
improving as the receiver approaches the station. Other than 


knowledge of the station location and the dot and dash 
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Figure 12. Polar Diagram Of Consol Pattern. (USNO, 1972, p. 
7-29) 
position correlations, little is required to use Consol. 


(USNO, 1972, p. 7-5) 


F. VOR/DME AND TACAN 

VOR/DME stands for VHF Omni-directional Ranging/ Direction 
Measuring Equipment and TACAN stands for Tactical Air 
Navigation. It is, in effect, a radio beacon system for 
aircraft along heavily travelled airline routes and between 


military airfields and ships. The principle of operation of 
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TACAN is the same as for VOR/DME except that TACAN uses higher 
frequencies and is more accurate (Beck, 1971, pp. 26-28). Some 
VOR stations are equipped to respond to TACAN DME frequencies 
and are called VORTAC stations. VOR/DME uses VHF frequencies 
between 108 MHz and 118 MHz and TACAN uses UHF frequencies 
between 960 MHZ and 1,215 MHZ (Logsdon, 1992, p. 36). 
Consequently, both systems are line of sight only. 

VOR employs two different signals working together in 
partnership: a narrow-beam rotating "lighthouse" transmission 
coupled with an "blinking" omni-directional pulse. As the 
lighthouse transmission sweeps past magnetic north, the 
blinking pulse is transmitted. Receiver equipment marks the 
time of receipt for the first pulse and then the time of 
receipt for the lighthouse sweep. Since the lighthouse sweep 
rotates at a fixed 30 rev/sec, it is a simple matter to 
compute the bearing from magnetic north. VOR provides only 
bearing information. (Logsdon, 1992, pp. 36-38) 

For range from the VOR station, the aircraft transmits a 
series of interrogation pulses to the DME at the VOR station 
which then transmits a reply a fixed time after receipt. 
Equipment on the aircraft computes the slant range by dividing 
the time from interrogation to the time of receipt of the DME 
response by two and multiplying by the speed of light. Figure 
13 illustrates this operation. (Logsdon, 1992, pp. 36-38) 

At the extreme range of receipt, VOR can be off by as much 


as three degrees. However, as the range to the transmitter 
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decreases, the error is reduced. Overall, VOR/DME TACAN 
positioning accuracies range from approximately 200 to 600 
feet (Logsdon, 1992, p. 36). The obvious advantage over a 
radio beacon system is lack of need for RDF equipment on board 


the aircraft. 


G. OMEGA 

The Omega navigation system was designed to provide 
worldwide coverage and to be accurate. It consists of only 
eight stations, each uSing a base frequency of 10.2 KHz, and 
transmitting at 10 Kw power. This makes it possible to receive 
at least three and usually four stations anywhere on earth 
with an accuracy of approximately one nautical mile. Omega’s 
principle of operation is similar to Decca in that it is a 
hyperbolic radio navigation system that uses the phase 
comparison of two CW transmissions to obtain a line of 
position. It is different in that any two stations can be used 
to obtain a line of position and with special receiving 
equipment, range from the two stations for a fix. (Bowditch, 
1984, pp. 1016-1036) 

Two 10.2 KHz continuous waves transmitted exactly in phase 
but travelling in opposite directions produce a series of 
hyperbolic lanes called Omega lanes. Each lane is eight 
nautical miles wide which corresponds to the half-wavelength 
of the frequency. Phase comparison determines the receiver 


position within the lane which is expressed in terms of 
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Figure 13. VOR/DME Principles Of Operation. (Logsdon, 1992, 


centicycles (cec) or centilanes (cel). These are defined to be 


0.01 of the width of the lane. Figure 14 shows two Omega lanes 


produced by the 10.2 KHz frequency. (Hobbs, 1981, pp. 259-262) 
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Figure 14. Two Half-Wavelength Omega Lanes Produced By The 
Phase Comparison Of The 10.2 KHz Signals. (Hobbs, 1981, p. 
260) 


Lane ambiguity is resolved by multiplexing the basic 
frequency with three others at 11.05 KHz, 11.33 KHz and 13.6 
KHZ over a ten second period according to the commutation 
pattern shown in Figure 15. Within the receiver, the 
difference frequencies of 0.283 KHz, 1.133 KHz and 3.4 KHz are 
extracted which are 1/36, 1/9 and 1/3, the frequency of the 
10.2 KHz signal, respectively. The 3.4 KHz signal is used to 
establish a broader phase lane that encompasses three Omega 
lanes. Likewise, the 1.133 KHz signal phase lane encompasses 
three of the 3.4 KHz phase lanes and so on. By performing 
successive phase comparisons as demonstrated in Figure 16, the 
operator can easily determine within which Omega lane his 


receiver is located. (Hobbs, 1981, pp. 259-262) 
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H. SHIP’S INERTIAL NAVIGATION SYSTEM (SINS) 

Although it requires no external input, SINS is still an 
extremely useful navigation system. It is found primarily on 
submarines and aircraft carriers. Inertial navigation is... 


"... the process of directing the movements of a vessel 


based on sensed accelerations in known spatial directions 
by means of instruments that mechanize the Newtonian laws 
of motion, integrating such accelerations to determine 
velocity and position." (Hobbs, 1981, p. 281) 


In essence, an inertial navigation system such as SINS is a 


highly accurate dead reckoning computer. 
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Figure 16. Three Successive Omega Phase Comparisons For Lane 
Resolution On Surface Ships. (Hobbs, 1981, p. 261) 


SINS mounts two accelerometers ona platform stabilized by 
a system of three gyros in such a way that it is constantly 
coincident with a plane tangential to the earth’s surface. The 
two accelerometers are continually oriented in a north-south 
and east-west direction; hence, they are sensitive only to 
horizontal north-south and east-west accelerations. 
Integrating the accelerations with respect to time gives the 
component velocities from which the ship’s true velocity may 
be computed. (Hobbs, 1981, p.281-282) 

Because of the cumulative affect of a number of possible 
sources for errors such as friction in the gyro supports and 
disturbances caused by the daily rotation of the earth, SINS 
accuracy degrades over time and must be updated from external 
sources. One improvement to SINS has been the Electro-Static 
Gyro (ESG) which is used to reset or update the SINS. Simply, 


the ESG is a one centimeter diameter beryllium sphere spinning 
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at 216,000 rpm in a near perfect vacuum. This rotor is 
supported solely by an electrostatic field. The sphere is thus 
freed from the classical gyro bearing friction as well as many 
of the associated random torques that a mechanical support can 
introduce. The ESG extends the interim period between external 
updates by a factor of six. (Kayton, 1990, pp. 194-206) 

Another improvement in SINS technology is the Ring Laser 
Gyro (RLG). It differs from the traditional concept of a gyro 
in that it contains no spinning mass. It is a closed geometric 
laser path (usually triangular such as the one in Figure 17) 
centered on an expected spin axis. Identically phased laser 
beams are continuously generated, which travel in opposite 
directions around the closed path. Any rotation about the spin 
axils causes an apparent phase difference in the two beams at 
a measurement point at one vertex of the triangle, 
proportional to the speed of rotation. When combined with the 
output of three mutually orthogonal accelerometers, position 
and velocity are easily computed. RLG’s increase the precision 
over an ESG monitored SINS by an order of magnitude. (Hobbs, 
1981, p.281-282) 

Finally, in an effort to further improve accuracy and 
reduce size and weight, the latest development is the Fiber 
Optic Gyro. Its principle of operation the same as RLG’s 
except the laser path can be tightly wound in a spool and 
mirrors are not required. The technology is very young and 


promising and is still being perfected. 
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Figure 17. Single Axis Ring-Laser Gyro. By Monitoring The 
Interference Pattern Caused By The Two Counter Directional 
Isophase Light Beams, Extremely Small Rotations Can Be 
Measured. (Logsdon, 1992, p. 107) 


I. SHIPBOARD DOPPLER SONAR SYSTEM 

Another navigation system that does not provide absolute 
position as much as it provides accurate inputs to a dead 
reckoning system is the Shipboard Doppler Sonar System. 
Because of its extreme accuracy in measuring very slow 
velocities, it is found on very large ships where ship’s speed 
over ground is important. For example, when mooring to a pier, 
the ship could crush the pier if the contact velocity were 
greater than a Knot (Hobbs, 1981, p. 285). 

It operates by measuring the doppler shift of two or more 


sonar beams in order to compute speed over ground. One such 
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system by Raytheon, called Janus, uses two fixed hull mounted 
transducers that transmit pulsed sonar signals in two beams, 
one forward and the other aft, at about a 30° angle from the 
vertical. Figure 18 illustrates a similar configuration that 
uses four transducers. The former of the two produces fore and 
aft speed accuracies of 0.1 kt, depth accuracies to the 
nearest foot, meter or fathom and distance accuracies to the 
nearest hundredth of a nautical mile. (Kayton, 1990, pp. 188- 


193) 


J. TRANSIT 

Transit Navigation Satellite System, sometimes called just 
Transit or SatNav, can trace its conceptual origins directly 
to the launch of Sputnik I. It consists of five satellites in 
polar orbits at an altitude of 580 nm with periods of 
revolution of 107 minutes as shown in Figure 19. The rest of 
the system includes ground tracking stations, a computing 
center, an injection station and receivers. Its geometry is 
configured such that every satellite comes within range of 
every position on earth at least twice a day at 12 hour 
intervals. The average time between fix opportunities 
approximately 100 minutes at the equator and a half an hour 
near the poles. The difference is due to the intrinsic nature 
of polar orbits. (Bowditch, 1984, pp. 1066-1095) 

Transit employs two frequencies in the UHF band of 150 MHz 


and 400 MHz. It operates by using an integrated doppler method 
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Figure 18. A Dual-Axis Four-Beam Janus Sonar Array Of The 
Raytheon Doppler Sonar System. (Hobbs, 1981, p. 285) 


where the slant range from the receiver to the satellite is 
computed by fitting the doppler phase shift curve of the 
satellite to look up table. Each satellite transmits a phase 
modulated message that repeats every two minutes containing 


the satellite time, ephemeris and other information from which 
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Figure 19. The Distinctive "Bird Cage" Polar Orbits Of 
Transit Satellite System. (Hobbs, 1981, p.275) 


the receiver calculates the receivers position in space and 
computes the longitude and latitude. Within in each pass of 
the satellite, the receiver needs to collect at least three 
successive messages for a fix. Four are preferred and seven 
are optimal. The number receivable is dictated by the height 
above the horizon achieved by the satellite. (Bowditch, 1984, 


pp. 1066-1095) 


eal 


The accuracy of Transit is affected by ships dead 
reckoning accuracy and speed of advance. For a stationary 
ship, it is common to achieve accuracies of 35 meters made 
possible by the shorter wavelength of the UHF signals 
(Bowditch,’ 1984, pp. 1067). Sit is Signifiteant to notemwtnam 
while this system is world wide and has improved positioning 
accuracies of better than an order of magnitude over systems 
discussed so far, it is comparable to Omega in its coverage 
and worse in its availability to provide positioning 


information. 


K. GLOBAL POSITIONING SYSTEM (GPS) 

With the June 1993 launch of the twenty-fourth satellite, 
NAVSTAR’s Global Positioning System (GPS) became fully 
operational in late 1993. GPS provides real time three 
dimensional positioning and timing information anywhere in the 
world. It consists of 24 uniformly spaced satellites in six 
orbital planes at an altitude of 20,180 Km as depicted in 
Figure 20. Like Transit, it uses two UHF frequencies at 1,575 
MHZ and 1227 MHz simultaneously upon which constellation 
ephemeris data, atmospheric propagation correction data and 
satellite clock error are phase modulated. (Wells, 1987, pp. 
4.00-4.11) 

The principle of operation for GPS is to calculate the 
range from four different satellites using the broadcast 


message, and to compute the longitude, latitude and altitude 
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Figure 20. GPS Constellation With 24 Satellites Ensures At 
Least Twelve In View At All Times Anywhere In The World. 
(Hobbs, 1981, p. 288) 


of the simultaneous solution of the resulting ranges. Because 
the system is operated by the Department of Defense for 
military use, the lower of the two frequencies is encrypted 
and altered to deny any potential enemy real time positioning 
accuracy for medium and long range targeting. St iis 
accuracies uSing just the higher frequency are approximately 
100 meters (Logsdon, 1992, p.64). Using methods to be 


explained in the next chapter, this can be improved by more 
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than an order of magnitude and in non-real time applications 
such aS surveying, mind boggling accuracies of a millimeter 
spherical error probable (SEP) can be achieved using phase 


comparison techniques (Logsdon, 1992, p.64). 


L. SUMMARY OF CONTEMPORARY ELECTRONIC NAVIGATION SYSTEMS 

It is easy to see why, after surveying the various 
electronic navigation systems available, GPS invites the most 
enthusiasm and excitement from all users. It is a quantum 
improvement over anything else available. There are other 
navigation systems not discussed here such as Microwave 
Landing Systems (MLS), the French Argos system and the Russian 
GPS imitation GLONASS. All suffer from one or more limitations 
in range, coverage, availability or accuracy. What is 
important in this review is an appreciation of the various 
contemporary electronic methods and their method of operation 
including their particular limitations. In the next chapter, 
the GPS principle of operation is more closely examined as 


well as the various methods for its use. 
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III. GLOBAL POSITIONING SYSTEM 


GPS operational and technical characteristics must be 
understood in detail in order to follow the reasoning for GPS 
receiver selection and the concept of employment in 
interfacing them with the CPU. In the following sections, a 
comprehenSive review of the User Segment of GPS is provided. 
It covers the salient points of signals, coding, message 
content, message structure, calculations, and various biases. 
The chapter ends by using this information in a description 
and discussion of Differential GPS. 

GPS is functionally divided into three segments: Control 
segment, Space segment and User segment. The Control segment 
exists to monitor, update and operate the whole system. The 
Space segment consists of the 24 satellites. These two 
critically important segments, however, are of no concern to 
the user and little 1s required to be known of them. The User 
segment is the most important, for understanding its operation 
will greatly influence the user’s choice of receiver and 


subsequent method of employment. 


A. USER SEGMENT RECEIVER HARDWARE DESCRIPTION 
Currently, receivers cost as little as a few hundred 
dollars to as much as $50,000 or more. Size and weight also 


vary greatly from receivers as small as a package of 
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cigarettes weighing approximately the same to larger than a 
car battery and weighing only slightly less. Shape and type of 
antenna can also impact the effectiveness of the receiver. 
Figure 21 illustrates the most common designs of antennae 
found on the market. 

One type of antenna called a microstrip is relatively thin 
(approximately 1/10" to 1/2" thick) and can be designed to 
receive one or both of the frequencies that GPS satellites 
transmit. This design is ideal for use in this project where 
low drag is desired on an aircraft. However, signal reception 
from satellites low on the horizon usually experience a much 
lower signal-to-noise ratio due to the lower gain of the 
antenna at those low angles caused by its thin design. 

Other factors that distinguish receivers from one another 
are the choice of position solution algorithm and interface 
communication protocols. All of the above characteristics are 
determined and limited by the state of the art of the 
technology and by the purpose for which it was designed. 

Receivers may be loosely grouped into four’ broad 
categories based on capability (Logsdon, 1992, p. 66): 

1. Number of channels and sequencing rate. 
2. Access to selective availability signals. 
3. Use of available performance enhancement techniques. 
4. Computer processing capabilities. 
The majority of user receiver sets sold are continuous 


tracking or parallel receivers. They typically have from four 
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Figure 21. Four Common Types Of GPS Antennae. (Wells, 1987, 
mee 7 .O2) 


to twelve channels. There are also single channel receivers 
capable of providing relatively crude positioning. To obtain 
data from multiple satellites, they employ one of three 
different sequencing speeds: slow, fast or multiplexed (very 


fast) as shown in Figure 22. (Logsdon, 1992, p. 66-67) 
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Figure 22. Multiplexing And Sequencing Channels For GPS 
Receivers. (Wells, 1987, p. 7.13) 


Selective Availability (SA) is the primary means by which 


the United States Department of Defense denies the potential 
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terrorist or enemy accurate use of GPS for targeting 
navigation for missiles. GPS messages are encoded using two 
codes. SA alters and encrypts the message content of the more 
accurate Precision (P) code carried on both frequencies. The 
less accurate Course and Acquisition (C/A) code message is 
called the Standard Positioning Service (SPS) and is available 
to anyone on the higher frequency only. The P-code information 
is called the Precise Positioning Service (PPS). Only those 
receivers with a decoding ip (tightly controlled by the 
DOD) can receive PPS information. (Logsdon, 1992, p. 67) 

SPS receivers can still achieve excellent accuracies by 
using performance enhancement techniques. One such technique 
capable of centimeter accuracy uses interferometry. This 
process needs very expensive antennae and receivers, and the 
data has to be post-processed. Surveyors are the most common 
users with this accuracy requirement and do not need the real 
time positioning information so denial of PPS is not a 
Significant handicap. Another capability that enhances SPS 
accuracy in near real-time uses a process called Differential 
GPS (DGPS). Simply, DGPS employs two receivers: one stationary 
at a known location and one mobile. The stationary receiver 
computes the difference between its known location and its GPS 
derived location. It then transmits those differences to the 
mobile receiver which uses them to correct its own GPS derived 
position. Since the Selective Availability errors are the same 


over baselines as large as 600 nm, accuracies approaching and 
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exceeding those of PPS receivers can be achieved. DGPS is 
explained more completely later. (Logsdon, 1992, p. 67) 

The last category of receivers is distinguished by 
computer processing capability. Mostly esoteric, this is where 
the manufacturer takes the GPS message and calculates the 
output message accounting for all the nuances of the solution 
algorithm. These include such items as whether or not to use 
just four satellites or more when they are available, whether 
or not to account for ionospheric and tropospheric delays 
(sometimes modeling these delays can introduce more errors 
than they eliminate), when to shift to other satellites that 
will reduce the geometrical error, etc. In addition, the 
receiver can be enhanced through the use of a Kalman Filter to 
help smooth the position and time data as well as reduce 


receiver processor time. (Logsdon, 1992, p. 68) 


B. GPS SIGNALS AND CODING 

GPS satellites continuously transmit two frequencies: Ll 
at 1575.42 Mhz and L2 at 1227.60 Mhz. The respective 
wavelengths for Ll and L2 are approximately 20 cm and 25 cm. 
Each of these is modulated with the GPS message at 50 bps. 
Since each satellite uses the same frequencies to send its own 
unique message, the receiver must be able to distinguish one 
Signal from the other. Furthermore, the transmitting power of 
each satellite so is low that by the time the signal reaches 


the receiver, its signal-to-noise ratio is less than the 
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ambient noise level. GPS overcomes these obstacles by using 
spread spectrum and code division-multiple access techniques. 
As mentioned earlier, the messages are encoded using two codes 
called C/A code and P-code. The chipping rate of the C/A code 
and P-code is 1.023 Mbps with a virtual wavelength of 290 m 
and 10.23 Mbps with ae virtual wavelength of 29 wm, 
respectively. Each satellite uses its own unique code with 
which the receiver may auto-correlate a duplicate and thus 
phase lock on to it. The combination of the transmission, 
coding and message chipping rates provides for a very robust 
Signal with which it is very difficult to interfere. (Wells, 
1987, pp. 4.5 - 4.10) 

The C/A code uses a_ 1023 bit Gold code with a period of 
one millisecond. The Gold code was chosen because of its 
length and because of the property that Gold codes are nearly 
orthogonal. It is also very easy to replicate and with its 
short period, easy to achieve correlation with itself. A GPS 
Signal can only be received if the code assignment is known 
ahead of time. The Gold code assigned to each satellite is 
published in the Interface Control Document ICD-GPS-200 
published by Rockwell Space Systems Division. (Wells, 1987, p. 
6.07) 

The P-code uses a much longer code with 
235,469,592,765,000 bits and has a period of approximately 267 
days. It is divided into 32 seven day segments, and a segment 


is assigned to each satellite. The DOD further encrypts the P- 
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code to produce the Y-code to deny its use to unauthorized PPS 
users. The encryption and length of the P-code make it 
extremely difficult to correlate to without keying information 
and decryption hardware. (Wells, 1987, p. 6.07) 

The GPS message iS phaSe modulated onto its respective 
codes which in turn are phase modulated onto their respective 
carrier frequencies. As mentioned earlier, the L1 frequency is 
modulated with both the C/A code and the P-code. This is 
accomplished by phase modulating its quadrature component with 
the P-code. This feature primarily allows the military to 
employ PPS redundancy in certain applications. The table in 
Figure 23 compares the important points of the code and data 


features. (Wells, 1987, pp. 6.3 - 6.8) 


C. GPS MESSAGE CONTENT 

A GPS receiver determines its position by computing the 
position ot four or more satellites within view and then 
measuring its range from each. By solving four simultaneous 
equations, the receiver calculates its coordinates in an Earth 
Centered Earth Fixed (ECEF) cartesian coordinate system. The 
system chosen is defined to be consistent with the World 
Geodetic Survey of 1984 (WGS-84) which places the origin at 
the intersection of the Z axis (through the earth’s center of 
mass and parallel to the north pole of the Conventional 
Terrestrial Pole (CTP)) and the X-axis (the intersection of 


the WGS-84 Reference Meridian Plane parallel to the zero 
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SPECIAL Easy to Slightly more accurate. | Provides handover from 
paen TIES acquire Resistent to jamming | C/A to P code, ephemeris 
and spoofing. | data and clock correction. 
Rejection of multipath. 











Figure 23. GPS Code Comparison. (Applied Research 
Laboratories, 1993, p. 2.2) 


meridian plane, and the plane of the CTP’s equator). The Y- 
axils completes the ECEF right hand orthogonal coordinate 
system measured in the equatorial plane 90° to the east as 
shown in Figure 24 (Applied Research Laboratories, 1993, p. 
2.1). Having defined the reference frame for GPS position, it 
is necessary to understand a few more concepts in describing 
position within an ECEF system. 

The GPS receiver translates the ECEF coordinates into 
geodetic coordinates for the user. The equations to translate 


coordinates from ECEF to geodetic coordinates and vice-versa 
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BIH - Defined CTP 
tWGS B4 


Earth's Center 
of Mass 


BIH - Defined 
Zero 
Meridian 





Figure 24. GPS Earth Centered Earth Fixed Coordinate System 
As Defined By The WGS 84. (BIH Stands For Bureau 
International De L’Heure) (Applied Research Laboratories, 
LOO epee oils) 

are provided in Appendix A. By using the reference ellipsoid 
model of the world defined by WGS 84, GPS is capable of 
relating pin point positions between continents. Before the 
World Geodetic Surveys of 1972 and 1984, inter-continent 
surveying was virtually impossible. Cartographers could only 


establish local datums to which all surveyed points were 


referenced. These datums were quite accurate for their 
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respective geographical regions at best. One of the 
accomplishments of WGS 84 was in relating all the different 
local datums to the ECEF coordinate system. By modelling the 
datum as a point of tangency of an ellipsoid as shown in 
Figure 25, the ellipsoid origin could easily be referenced to 
the origin of the ECEF coordinate system. With the advent of 
WGS 84 and GPS, the importance of the local datums has been 
diminishing. Most GPS receivers contain a library of the major 
local ellipsoid references to the WGS ECEF coordinate system 
so charts and maps that use different local datums can be 
employed with GPS. 

Once the reference frame is established and understood, 
satellite position can be described using the six classical 
Keplerian orbital elements. They are briefly described in 
Figure 26 but they are not sufficient for the most accurate 
description of the satellite orbit. This is because the 
predicted position of a satellite is the result of least- 
squares curve fitting which is not valid for an entire orbit. 
Thus, it is necessary to include several other parameters to 
more accurately model satellite positions in space and time. 
All the parameters are listed in Appendix B. All told, there 


are 17 specific parameters included in the GPS message. 


D. GPS CALCULATIONS 
The equations required to determine the ECEF coordinates 


of the satellites from the ephemeris are presented in Appendix 
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Local ellipsoid 





Figure 25. Two Dimensional Illustration Of A Theoretical 
Ellipsoid Tangent To The Geoid At An Arbitrary Datum P. 
(Leick, 1990, p. 189) 


B. If the receivers and satellites could be perfectly 
synchronized with GPS time, three parameters would be unknown: 
X, Y, Z. It would then be possible to determine position using 
only three satellites. However, all receivers use inexpensive, 
lightweight and comparatively inaccurate quartz clocks and 
each has a bias from GPS time for which there must be an 


accounting. Therefore, for three dimensional position, there 
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THE SIX KEPLERIAN ORBITAL ELEMENTS 


THE SIX KEPLERIAN ORBITAL ELEMENTS 
a@lQawt 


a@€: SEM! MAJOR AXIS @: ORBITAL ECCENTRICITY 


THE “OBLATENESS” 
THE HALF - LENGTH OF THE ORBIT 
| | OF THE ELLIPSE a ra- 


~4—-——. 78 ——_ —» fate 
fp+fs 


1: ORBITAL INCLINATION Q: ASCENDING NODE 


THE ANGLE BETWEEN 
THE EQUATORIAL PLANE 
AND THE ORBITAL 
PLANE 


@®: ARGUMENT OF PERIGEE 
THE LOCATION OF THE THE TIME OF PASSAGE 
PERIGEE POINT OF THE POINT OF 
CLOSEST APPROACH 


Figure 26. Six Classical Keplarian Orbital Parameters. 
(Beagqsdon, 1992, p. 143) 


THE LONGITUDE OF THE 
ASCENDING NODE 
(EQUATORIAL CROSSING) 





are four unknowns and four simultaneous equations as shown in 


Equation (2): 


Pe Se a) ee i Ve) a ee CUAL. 
Pee, (ce xe) eniy- t-ny ene tae oe Ae. Sy 
ek Jw) Oca, eer oe ms a eS Toons 
eh oy cai aan ee (RS Og I es a 

where P,” is the pseudorange for the «th epoch in time, # is 


the satellite number, c is the speed of light and At,is the 


receiver clock bias from GPS time. 
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Parenthetically, it should be noted that although GPS time 
is synchronized within one microsecond of UTC, it does not 
recognize UTC leap seconds. The last simultaneous epoch of GPS 
and UTC time was in January 1980. Consequently, GPS time 
differs from UTC by the integer number of leap seconds since 


then. 


E. GPS MESSAGE STRUCTURE 

The entire GPS message format is illustrated in Figures 27 
through 29. It is defined by ICD-GPS-200 written by Rockwell 
International who was the prime contractor for GPS. The basic 
message unit 1s called a Frame. It is subdivided into five 
subframes which are in turn subdivided into ten words each. A 
word has 30 bits. For positioning, a receiver requires only 
subframes one through three. Those contain the necessary 
ephemeris and clock corrections. However, the designers saw a 
need to provide the user information about all the satellites 
in the entire GPS constellation. They used subframes four and 
five to do this. Since this "Almanac" of extra information 
exceeded the size of those two subframes alone, it is 
transmitted using subframes four and five of 25 successive 
message frames. Each pair of subframes four and five of the 25 
successive frames is referred to as a page. The entire block 
of 25 pages and subframes one through three is considered a 


Master Frame. Figures 27 and 28 illustrate this best. At a 
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data rate of 50 bps, a frame needs 30 seconds and a master 


frame requires 12.5 minutes for complete transmission. 


BASIC MESSAGE UNIT IS ONE FRAME (1500 BITS LONG) 


+—______ 30 sec 


1 FRAME = 5 SUBFRAMES ee caie 


NN 
1 SUBFRAME = 10 WORDS [1/2 [3]4]5]6/7 1/8 [9 [10 


1 WORD = 30 BITS SOC 
0.02 sec 


ONE!I MASTER . Subframes 4 and § 
S FRAME |INCLUDES te aa ene 


ALL 25 PAGES OF SUBFRAMES 4 & § 
= 37,500 BITS TAKING 12.5 MINUTES 





Figure 27. GPS Message Format Breakdown. (Rockwell, 1987, p. 
vets ) 


F. GPS BIAS 

There are numerous factors that affect GPS performance 
accuracy. Selective Availability is primary among these. While 
the exact nature of Selective Availability is classified, it 
is sufficient to know that it denies single receiver position 


accuracy of better than 100 meters 95% of the time. In 
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addition, the Ionosphere and Troposphere affect the time of 
receipt by slowing down the signal. If the receiver could 
receive two frequencies, these atmospheric errors could be 
determined and eliminated. For the receivers that are capable 
of only single frequency reception (most receivers), 
compensation for these errors iS achieved by using 
mathematical models which adequately reduce their affects for 


most positioning purposes. 


The GPS Data Message Content 


yma cies 


Flags (L2 code & data; week #; , : 
satellite accuracy and health) Alinanac for satcllites 1-24 (pages 1-24) 


Age of data Health of satellites 1-24 (page 25) 


Satellite clock correction coefficients 


Almanac for satellites 25-32 (pages 2, 3, 4, 5, 7, 8, 9, 10) 
Ionospheric model, and UTC data (page 18) 

Antispoof flag - 32 satellites (page 25) 

Satellite configuration - 32 satellites (page 25) 

Health of satellites 25-32 (page 25) 





Figure 28. GPS Message Content Organization Within The 
Message Format. (Applied Research Laboratories, 1993, p. 
2) 


Signal reflections, commonly referred to as multipath, are 


another source of position errors. Multipath interferes with 
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the fundamental concept of measuring the time from 
transmission to the time of receipt in order to determine 
range from the satellites. Fortunately, with adequate 
filtering and by paying careful attention in choosing the 
antenna, almost all the effects of multipath may be 


eliminated. 


The GPS Data Message 
50 Bits Per Second Data Rate 


Bit No. > 0 30 60) 
Subframe Telemetry Handover 


Words Words Clock Correction 


300 330 60 
Subframe Telemetry Handover 


Words Words Ephemenis 


630 660 


600 
Subframe Telemetry Handover 
Words Words 


Ephenicris 


900 
Subframe | Telemetry Handover 


Words) Words (Multipiex) Message (Changes Through 25 Frames) 


1200 4230 1260) 1Sin} 


Subframe Telemetry Handover AjJmanac/Healih Status 30) Sce** 
E C 


* Format of Frame Twenty-Five Changes. 
°° 32.5 Minutes Before The Entire Message Repeats. 





Figure 29 GPS Message Organization. (Applied Research 
Laboratories, 1993, p. 2.2) 


Finally, satellite geometry relative to the receiver plays 
a role second only to Selective Availability in determining 
the precision achieved. The measure of "goodness" of satellite 
geometry is called Geometric Dilution Of Precision or GDOP. 


Optimal GDOP (the lower the GDOP number the better) is when 


oypk 


three of the satellites are on the horizon equally spaced 120° 
apart and the fourth is directly over head as shown in Figure 
30. Since the constellation is designed such that there are 
usually 12 satellites in view at all times, a low GDOP of two 


or three can regularly be achieved. 
o = DOP :s6 


Positioning Geometry Measurement 
accuracy (Dilution of accuracy 
Precision) 


aH &§ 


POOR GDOP GOOD GDOP 
satellites bunched (ideal case) 
together ¢ one satellite overhead 
¢ 3 on horizon, 
120° apart in azimuth 





Figure 30. Geometric Relation Between GPS Satellites Affects 
The Position Accuracy. (Wells, 1987, p. 4.22) 


G. DIFFERENTIAL GPS 
SPS does not guarantee accuracies less than 300 m. It is 


possible to achieve accuracies two orders of magnitude better 


ore 


through a process called Differential GPS but it requires at 
least two receivers. The concept is simple and generates the 
greatest amount of promise in commercial applications of GPS. 
By placing one of the two receivers in a pre-surveyed location 
where its geodetic coordinates are known as precisely as 
possible, it can compare its true position to the position GPS 
provides and calculate the differences between latitude, 
longitude and altitude. If the resulting error signs are 
changed, the errors become corrections that may added to the 
second receiver’s GPS position for a greatly improved position 
accuracy. Over baselines of up to 600 nm, the errors are 
approximately the same. 

In practice, there are a few restrictions in employing 
DGPS. The most important is that the position results from 
both receivers must be derived using the same satellites or 
else the improvements in position are diminished by the 
introduction of new errors. The other restriction is that the 
corrections should be applied at the earliest point in the 
algorithm used to determine position coordinates in order to 
Minimize the propagation of errors. The reasons for these 
restrictions will become clear after a brief description of 
the algorithm used to compute DGPS corrections. 

The stationary GPS receiver is receiver is required to 
send corrections to the mobile receiver. It would be simple to 
subtract the true longitude, latitude and altitude from the 


geodetic position obtained from the GPS satellites. However, 
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the errors introduced from SA, satellite geometry and 
atmospheric interference affect the psuedoranges from each 
satellite in varying degrees. One satellite psuedorange will 
always have a greater error than the next. Since the GPS 
position is calculated using a least squares reduction of the 
four simultaneous equations listed in equation two, all the 
individual psuedorange errors will be "blended" thereby 
permitting the greatest error to act on all the psuedoranges 
and the final position solution. 

In order to avoid this problem, the stationary receiver 
computes the expected position of and range from each 
satellite (based on the time of travel of the individual 
Signal) and compares it to the position and range that the GPS 
message asserts in its ephemeris. The difference between the 
two becomes a psuedorange correction. At least four of these 
corrections are transmitted to the mobile receiver and applied 
to the corresponding psuedoranges obtained by that receiver. 
Then the new psuedoranges are used in equation two to 
calculate the longitude, latitude and altitude of the mobile 
receiver. Since the individual corrections are applied before 
equation two, the final position is more accurate. 

In practice, GPS makes all of its calculations using the 
Cartesian Coordinate ECEF system so the stationary receiver 
must translate its true position from geodetic coordinates to 
ECEF coordinates. In a similar manner, the mobile receiver 


must translate its corrected position from ECEF coordinates to 
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geodetic coordinates. To summarize, the following steps are 
involved in improving the GPS position using DGPS: 


1. Convert the true geodetic position of the stationary GPS 
receiver to ECEF coordinates. 


2. Compute the satellite positions from the ephemeris. 

3. Compute the expected satellite positions by multiplying 
the time from transmission to time of receipt by the speed 
Oem velit. . 

4. Compare the received pseudoranges to the expected 
pseudoranges and calculate the pseudorange corrections for 


each satellite. 


5. Transmit the pseudorange corrections to the mobile GPS 
receiver and add to the received pseudoranges. 


6. Compute the mobile receiver’s position using the 
corrected pseudoranges. 


7. Convert the mobile receiver’s position from ECEF 
coordinates to geodetic coordinates. 
H. SUMMARY 
The User segment of GPS draws together the results of 
several decades of study and experimentation. Hardware design 
and miniaturization are pushing the state of the art. Signal 
construction and frequency selection optimize its immunity to 
noise and interference. Orbital kinematics referenced to an 
ECEF coordinate system simplify and standardize the 
information contained in the GPS mes:age structure. Even so, 
natural and man-made biases corrupt the integrity and accuracy 
of GPS. Through DGPS users can overcome the degradations 
caused by these biases and achieve real time accuracies 


unheard of heretofore. Understanding these principles upon 


oe. 


which the User segment is founded will enable the user to 
better select a receiver that suits his purpose and integrate 


1t into whatever application he so desires. 
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IV. GPS RECEIVER SELECTION AND INTERFACE WITH CPU 


A. MISSION REQUIREMENTS FOR GPS 

Before any serious search for a GPS receiver can be made, 
it is necessary to define the extent to which it will be used. 
These criteria come from the mission statement. This project’s 
goals are to integrate an Inertial Measurement Unit (IMU) with 
GPS and other non-INS sensors through a smoothing and control 
filter to achieve autonomous flight of an Unmanned Aerial 
Vehicle using off the shelf hardware and technology. As a 
result, several requirements are defined for the GPS 
receivers. 

To begin with, the position output should be frequent 
enough to provide regular corrections to the IMU. Also, the 
GPS receiver must actually update its position solution as 
often as possible to ensure accuracy in the information of its 
outputs. The receiver that is to be in the aircraft should be 
small and lightweight. It should be easily interfaced to the 
CPU, and it should provide the necessary information to 
conduct phase measurements of the incoming GPS signal. (This 
last feature should provide the positioning accuracy to auto- 
land the aircraft.) Finally, it must be readily available on 
the open market. Intrinsic in the mission statement is that it 


must be affordable given the project budget. With these 
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guidelines, an informed selection can be made in choosing the 


best GPS receivers for the job. 


B. GPS RECEIVER SELECTION 

Having reviewed the principles of GPS operation and 
variations in receiver design features, an extensive survey of 
216 different GPS receivers from 54 companies was conducted. 
In the January 1993 issue of GPS World magazine, many of the 
critical facets of each receiver were listed for comparison 
(Chan, 1993, pp. 52-65). In order to stay within budget wenwi 
those receivers less than $5,000.00 were considered. (The 
highest priced receiver was $55,000.00 with many between 
$10,000.00 to $20,000.00.) Among the remaining units, those 
that were capable of tracking six or more satellites 
Simultaneously were selected for comparison. This left the 
list at seventeen receivers. 

Of the remaining seventeen receivers, all were capable of 
the same approximate accuracy and information output rate at 
once per second. The one notable exception was Navstar’s XR5 
with an output rate at 1/10 second, but ultimately, it was not 
chosen because, at $1000.00, it cost three times as much as 
the Motorola PVT-6. Since all the receiver units performed at 
virtually the same level, the deciding factors came down to 
cost, weight, size and power consumption. Table I shows the 


top eight receivers with these features listed. 
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Table I. COMPARISON OF TOP EIGHT CANDIDATES FOR THIS 
MeO eC lS GPS RECEIVERS. »(CHAN,, 1993, PP. 52-65) 


GPS Weight Power 
Maker (Watts) 


Furuno 
Magellan 
Magnavox 
Motorola PVT-6 


XR5 


Novatel 


Trimble 





The Trimble price quoted includes a starter evaluation 
kit. The Motorola receiver, when all options were installed, 
was actually more expensive. By the criteria used, Trimble had 
the better unit would have been chosen over Motorola but they 


were very slow in answering inquiries for characteristics and 
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pricing information. When the decision was made, they were not 
considered. However, two factors make the choice of Motorola’s 
PVT-6 the best selection. One, they also make a portable 
differential GPS base station called a LGT-1000. It is 
conceivable, in the future, this project will choose to 
experiment with the specialized equipment. And two, Trimble 
has since stated they do not intend to pursue the small, 
Single receiver user market. Instead, they will specialize in 
industrial scale receivers for use at airports and in 


surveying. 


C. GPS RECEIVER DESCRIPTION 
1. Hardware 

The Motorola PVT-6 iS an original equipment for 
manufacturer (OEM) module. Although evaluation software is 
provided ina starter kit, it is not useful when the CPU is 
multi-tasking with several other processes. PVT-6 stands for 
Position, Velocity and Time with six channels. Figure 31 shows 
the dimensions of the unit without the antenna. The 
specifications are listed in Appendix C. The PVT-6 is capable 
of tracking six satellites simultaneously. It uses only the 
C/A coded message from the Ll carrier. The code tracking is 
carrier aided. It can be powered by a 12 Vdc unregulated or 5 
Vdc regulated source. (Motorola, 1992, p. 3-2) 

The PVT-6 is supplied with a low profile antenna of a 


microstrip design called a strip patch antenna. The antenna 
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WEIGHT: 4 OUNCES or 110 GRAMS 


& MOUNTING SCREW 
THREAD SIZE M4X0.8 


i (4 REQUIRED) 





Figure 31. Physical Dimensions Of The PVT-6.(Motorola, 1992, 
pe 3-3) 


gain pattern is provided in Appendix D. The Ll frequency is 
collected by the antenna, passed through a narrow band filter, 
amplified and sent to the receiver unit via an RG-58 coaxial 
cable. It provides the antenna 5 Vdc @ 25 mA from the receiver 
and returns the amplified Ll signal from the antenna. It 
should be noted the power loss along the cable should not 
exceed 6 dB at a frequency of 1575.42 MHz which equates to a 
maximum length of six meters for RG-58 coaxial cable. RG-58 is 


the recommended cable because it allows a sufficient length of 
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cable while achieving lighter weight due to the small gauge 
without the corresponding losses due to increased resistance 
that accompany even smaller gauge coaxial cables at that 
frequency. Antenna placement should position the patch plane 
to be as level as possible and have an unobstructed view of 
the horizon. This ensures a direct line of sight to all the 
visible satellites. Table II and Figure 32 provide the power 


and data connector pin assignments. (Motorola, 1992, p. D-1) 


Table II SERIAL PORT AND POWER PIN ASSIGNMENTS. (MOTOROLA, 
meo2., eric: 352) 


aie Sn nnn 
| Pin # | Signal Name| 


(Optional) +5 Vdc regulated or +12 Vdc 
+12V/+5V BATT unregulated for running Real Time Clock and 
retention of satellite ephemeris information 
stored in Static RAM memory. 
(Optional) +5 Vde regulated for power 
iz +5V MAIN requirements of the entire GPS receiver(+1l2 Vdc 
Signal not used). 










ee | RTN |Power supply (+5V or +12V) return. 
es ley memory (EPROM) programming voltage. 


+12V MAIN +12vde unregulated for power requirements of the 
Cnertre, GPsereceiver. 

1 PPS? (Option A) 1 pulse per second output. 

1PPS RTN (Option A) 1 pulse per second return. 








RS232 TXD Serial RS232 data output. 
I 
Lag RS232 RXD Serial RS232 data input. 


i 


10 /|\RS232 RIN 





|Signal return for RS232 signals. 





'Not used (used only by Motorola for updating of software). 
“1 Hz pulse train with rising edge coincident with UTC/GPS 1 second 
tices 


As Figure 33 shows, the Ll analog signal enters the 


receiver and is down converted to an intermediate frequency 
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ANTENNA DATA/POWER 
CONNECTOR __—_ CONNECTOR 





Figure 32. Side View Of The PVT-6 Showing The Pin Numbering 
Geaer-. (Motorola, 1992, p. 3-2) 


(IF). The IF is passed to the 6-channel code and carrier 
correlator section where it is digitized and split for input 
into six separate channels for correlation, filtering, carrier 
tracking, code tracking and signal detection. The output is 
synchronously routed to the MPU (PVT-6 CPU) where the 
information is processed and formatted. Also within the MPU is 
the RS-232 interface for full-duplex, asynchronous data 
communication with the host equipment. (Motorola, 1992, p. 2- 
2) 

The PVT-6 is configured to operate at 9600 bh-'1d with 
eight bit words, one start and stop bit and no parity. Its 
communications port is considered Data Communications 
Equipment with respect to the RS-232 serial communications 
standard. However, close examination of the available pins 
for serial communications reveals the presence of only the TX 


(pin 8), Rx (pin 9) and Gnd (pin 10) pins. This dictates there 
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_ DIGITAL SIGNAL PROCESSING 
ANTENNA 
CONNECTOR 


RF DOWN CONVERTER oe 


RS-232 DATA 


7 NON-VOLATILE 
MEMORY 


DATA/POWER CONNECTOR 





Figure 33. Functional Block Diagram Of The PVT-6. (Motorola, 
1992, p. 2-1) 


will be no modem line configuration (DTR and DSR) and at least 
no hardware flow control (CTS and RTS). In fact, there is no 
software flow control (XON and XOFF) either implying the Data 
Terminal Equipment (DTE- host equipment) must monitor the port 
continuously or miss valuable information. Likewise, the PVT-6 
must do the same. In practice, it receives all serial input 
into a two kilobyte buffer from which it reads 750 bytes per 
second. The buffer is large enough that it never overflows. 


(MoGOoErola 61992 pp ee) 
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Finally, the PVT-6 is provided with static random 
access memory (SRAM) for retention of the satellite ephemeris 
data. To prevent the loss of this information, an external 
power source such as a battery may be hooked up to the 
+12V/+5V BATT connection (pin 1) when the GPS unit is powered 
off. In addition, nonvolatile EEPROM is used for storage of 
custom operating parameters, almanac information and other 
information such as the receiver serial number and option 
list. (Motorola, 1992, p. 2-2) 

2. Motorola Binary Message Format 

The PVT-6 receiver is capable of providing its 
information in any one of three formats: Loran, NMEA-0183 and 
Motorola Proprietary Binary (Motorola, 1992, pp. 4-2 - 4-6). 
Both Loran and NMEA (National Marine Electronics Association) 
would severely restrict this project’s flexibility and 
accuracy, and so are not used. Motorola Binary format is 
logically organized and provides tremendous flexibility to the 
user to apply GPS to what ever purpose. 

Motorola’s Binary format has 47 different messages it 
sends and receives ranging in length from 7 to 294 bytes. Each 
message begins with "@@" (double ASCII 0x40) which is followed 
by a two letter, case sensitive message ID. Then, O to 287 
bytes of data follow, a byte size check sum and an ASCII 
carriage return-line feed (0OxOD Ox0A). The check sum is the 


exclusive-or of all the preceding bytes in the message. The 
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data that follows the message ID can be thought of as scaled 
floating point data stored in integer format. This avoids the 
floating point storage convention as well as saves space (only 
one byte is needed instead of four or more). The user must be 
very careful to observe the units of measurement provided in 
the Motorola GPS Technical Reference Manual. (Motorola, 1992, 
p= 2) 

For every command to the GPS receiver, there is a 
response usually using the same message ID. The receiver has 
no way to determine if the information it has received from 
the CPU is correct other than by verifying the several 
intrinsic properties of the message itself: 

e The message begins with "@@"., 
e The message ends with CR-LF. 


e The length of the message matches that expected for the 
message ID. 


e The checksum is valid. 


e The message data is within the expected maximum and 
minimum ranges. 


If the command input to the PVT-6 changing a parameter is 
invalid for any of the above reasons, it responds as expected, 
but the data values will not reflect the ordered change. It is 
prudent, therefore, to compare every input message with the 
subsequent PVT-6 response. (Motorola, 1992, pp. 4-2 - 4-3) 
Input commands may be of the type that change a 
particular configuration parameter of the receiver. Any 


command with a message ID that begins with a upper case "A" 
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falls into this category. Input commands may also be of the 
type that enable or disable the output of data or status 
messages. They are the commands that the CPU will use for 
measuring position, velocity, time, pseudorange and satellite 
ephemeris data. These output status messages include 12 that 
are the work horses of the PVT-6. Table III shows only nine 
because three of the commands are the reciprocals used by the 
CPU to talk to the PVT-6. The nine shown are listed in order 
of precedence for transmission from the GPS receiver to the 
CPU. (Motorola, 1992, p. 4-3) 

Table III DATA MESSAGE OUTPUT RATES SHOWN IN ORDER OF 


PRECEDENCE OF TRANSMISSION .20M THE PVT-6 TO THE CPU. 
(MOTOROLA, 1992, P. 4-4) 




















jj) DUTPUT. MESSAGE TYPE 224}, CONTINUOUS. (@=1.-255)/1 | TONE TIME (m-0) 2 = 


Visible Satellite Status When Vis data changes When requested 







es 


3. PVT-6 Operation 
The PVT-6 has some limitations that, if understood, 
may be overcome. Data is output once per second from the 


receiver. This information is calculated from the ephemeris 


6/ 


and time measurements of the previous second time tic called 
the measurement epoch. In order to provide the most accurate 
position given the latency of data, the position, velocity and 
time represent the best estimate of the data advanced one 
second. Figure 34 graphically illustrates this process. The 
first byte of information for a particular measurement epoch 
is transmitted anywhere from zero to 50 ms after the 
subsequent measurement epoch. Only the time for the Satellite 
Range Data message and the Pseudorange Correction message is 
not propagated forward. Careful consideration addressed to 
this latency issue in the design of any Kalman Filter or 
smoothing filter that integrates GPS position with the IMU and 
other sensors will improve performance. (Motorola, 1992, pp. 
4-13 - 4-17) 

To assist in synchronization of the GPS data with 
other sensor data, a one pulse per second (1PPS) signal is 
output using pins six and seven of the serial port. The pulse 
is approximately 200 milliseconds long with its leading edge 
corresponding precisely to the measurement epoch. The PVT-6 
has the option of additionally offsetting the 1PPS up to 
0.999999999 seconds. If there is a requirement for the 1PPS to 
be precisely synchronized with the UTC or GPS measurement 
epoch, the user may also offset the receiver measurement epoch 
as well as account for antenna cable delay. For this project, 


these latter conveniences will not be required. This is shown 
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Figure 34. Position/Channel/Status Message Data Output 
Latency. (Motorola, 1992, p.4-15) 


in the timing diagram in Figure 35. (Motorola, 1992, pp. 4-16 
~ 4-17) 

Finally, the PVT-6 suffers from periodic outages like 
most other receivers regardless of quality. These outages may 
be caused by temporary loss of satellite track do to obscuring 
by a building, tree or wing tip if the aircraft were in a 
banking turn. Sudden accelerations may also cause loss of 
track. In this event, the time for the receiver to reacquire 


the track and output data is important. As shown in Table IV, 
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Figure 35 Timing Diagram For The Epoch, 1PPS And Serial 
Output Data. (Motorola, 1992, p. 4-17) 
this time may take as long as 15 seconds. If the aircraft were 


attempting to land using GPS alone, this could be fatal. 


(Motorola, 1992, p. 4-7) 


Table IV NOMINAL REACQUISITION TIMES FOR THE PVT-6. 
(MOTOROLA, 1992, P. 4-7) 
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D. INTERFACING THE GPS RECEIVERS 


The two GPS receivers to be used for DGPS will be linked 
to the ground CPU and the UAV CPU through RS-232 serial ports. 
The big picture of the project showing where GPS fits in is 
shown in Figure 36. The ground control station will be 
interfaced through a 33 MHz Intel 80486DX personal computer. 
It will be the primary data flow conduit for the uplink and 
downlink as well as provide formatting and mass storage. The 
additional link in the ground station shown in Figure 36 is 
for redundancy in the event the main link fails and emergency 


control of the UAV is required. 
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Figure 36 Archytas Project Control/Communications Diagram 
Showing Where GPS Fits Into The Big Picture. 


al 


The UAV CPU uses off the shelf hardware that economizes on 
size, weight and price. This is done by using six computer 
cards and a passive backplane. The CPU is a 33 MHz Intel 
80486DX with 256 Kb double cache memory. Non-volatile memory 
for the operating system and programs 1S provided by a 2.88 Mb 
RAM/ROM Disk card instead of a relatively heavy hard drive. In 
addition, servo control and sensing are managed using an 
Analog-to-Digital (A/D) /Digital-to-Analog (D/A) card and a ten 
channel Timer/Counter card. Serial I/O is handled through the 
two standard serial communication ports on the CPU card (COM1 
and COM2) and an intelligent 8 port RS-232 serial I/O card. 
The GPS receiver communicates to the CPU through port 1 on 
this card (not to be confused with COM1 on the CPU card). In 
order save space, all the cards are connected to a common bus 
using a passive backplane. A simplified drawing of the final 
arrangement is shown in Figure 37. 

1. PCL=-744 8-Port Intelligent RS-232 Interface Card 

The UAV CPU needs multiple serial ports for other data 
inputs from the link modem, the IMU and the non-INS sensors. 
These multiple serial inputs are provided through Advantech’s 
PCL-744 8-Port Intelligent RS-232 Interface Card. 

The PCL-744 is an intelligent card in that it has its 
own processor and 64 Kb RAM buffer. It requires only one 
interrupt number for eight ports although each may configured 


independently. This configuration includes determining the 
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Figure 37 Simplified Illustration Of UAV CPU Hardware. 


size of the transmit and receive ring buffers and what type of 
flow control is used. These facts enable the UAV CPU to 
perform multi-tasking by allowing it to access the ports and 
process data without losing any information from the other 
ports arriving simultaneously. The PCL-744 takes care of that 
by multiplexing between its eight ports and placing the 
incoming data in buffers that can be accessed by the UAV CPU. 

The PCL-744 comes with a library of 40 user functions 
written in several high level languages including Microsoft C 
and Turbo C by Borland. This project uses the Turbo C routines 
for all serial communications. These functions allow the user 
to use the card in anyway desired for serial I/O. It is also 
necessary to ensure the TSR drivers are loaded prior to PCL- 


744 operation. They are called STD-DRV.EXE and 744-DRV.EXE and 
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should reside in the AUTOEXEC.BAT file to ensure their 
loading. 

The description of the software drivers for the GPS 
receivers that follows is the same in functional design for 
both the ground and UAV units, but there are _ slight 
differences in coding. This is due to the fact that the UAV 
receiver must communicate through the PCL-744, and the ground 
receiver communicates directly with the CPU through a standard 
DOS COM port. The port interface for the ground CPU uses a 
standard C library interrupt routine, and the port interface 
for the UAV that communicates through the PCL-744 uses the PC- 
ComLIB serial communication library routines provided with the 
Cala: 

2. GPS Software Driver 

In communicating with the PVT-6, all transmissions are 
without any flow control. It is necessary, therefore, to 
provide a storage buffer where the contents of the 
transmission may be stored for holding until the CPU is ready 
to handle the data. The information sent to the PVT-6 is put 
into a two kilobyte buffer in the receiver itself. All message 
handling 1S done by the receiver including parsing, 
verification and interpretation, and is transparent to the 
user. The user’s only responsibility is to ensure that the 
message sent to the receiver is properly formatted. However, 


when receiving messages from the PVT-6, parsing, verification 
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and interpretation must be done by the user. This allows 
flexibility for the user and is one of the minor reasons for 
selecting a Motorola GPS receiver. 

This section describes the software driver required to 
accept information from the GPS receivers. The software in 
Appendix E is written in modular format to allow flexibility 
and ease of handling. It is structured so it can be expanded 
and built upon. The following paragraphs describe the flow 
chart shown in Figure 38 which illustrate the concept of the 
GPS software driver in Appendix E and the structure of the 
software driver. 

The PCL-744 requires the eight ports to be initially 
configured using the vendors software. This configuration is 
stored in the PCL-744 driver and the card is thus reconfigured 
upon startup every time power is applied. This does not make 
the card ready for operation, though. Each port must be 
individually opened and set to the communication parameters 
required by its user. For instance, the PVT-6 uses 9600 BAUD 
Witeheeatn elght bit word, no parity and one start/stop bit. In 
addition, RS-232 modem control (1.e., DTR and RTS), hardware 
control (i.e., DTS and RTS) and software flow control (i.e., 
XON/XOFF) are set to off. This is all accomplished using the 
function "Open ports". At the conclusion of using the ports 
they should be formally closed and this 1s accomplished by 


using the function "Close ports". 
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Figure 38 Functional Diagram Of The GPS Software Drivers For 
Receiving Information From The PVT-6. 


The main functions that handle incoming information 
are called "Master gps" for the ground CPU and "Slave Gps "eae= 
the UAV CPU. The "Master gps" function first gets the buffer 
length of the message waiting in the Rx buffer for port number 
one on the PCL-744. It then creates a buffer of exactly that 
size and extracts the information from the Rx buffer and puts 
it into the temporary buffer just created. Next, it creates 
the input buffer which will contain the temporary buffer 
appended to the excess buffer contents. The excess buffer has 
the partial message left over from the previous call to 


"Master gps". With the input buffer created and the contents 
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of the excess and temporary buffers copied into it,the excess 
and temporary buffer memory is freed. 

The pointer to the input buffer is then passed to the 
function "Parse inbuff" where the first message is parsed from 
the input buffer and returned. If a complete message is not 
returned, the function places the partial message in the 
static excess buffer, frees the memory for the input buffer 
and exits the function "Master gps". If a complete message is 
returned, it is passed to the function "Msg verification" 
where it is checked for the correct length, check sum, ASCII 
prefix "@@" and hexadecimal suffix x0D x0A (ASCII carriage 
return and line feed). When it has been verified, the message 
is placed ina buffer to be passed out to the calling routine 
for transmission to the UAV GPS if it is a pseudorange 
correction. Otherwise, it 1s written to disk for post- 
processing. 

In the UAV GPS driver, the "Slave gps" function 
operates exactly the same way except after a message has been 
verified, it is put into a structure and scaled according to 
the format in the Motorola protocol for subsequent use by the 
Kalman filter. One of the more complicated structures for the 
ephemeris message has been written as an example by Motorola 
and incorporated in the header file GPSCONV.H. All message 
output from the UAV GPS is passed out to the calling routine 
to be down linked to the ground in binary format (to reduce 


downlink time) for post-processing. 
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Two other functions are mentioned but not yet 
completely written which play an important role. They are the 
"Initialize master gps" and "Initialize slave gps" functions. 
They will pre-configure the receivers for such things as true 
geodetic coordinates of the ground antenna, time, date, output 
rates, desired files to be collected etc. They should be made 
interactive to allow the user to change and set any 
parameters. The code framework for these functions is in the 
file GPSPORTS.H in Appendix E. 

Finally, the source file GPSAN.C with GPSCRNCH.H was 
used to post-process data collected for this thesis. It was 
written to take the Position/Channel/Status message from an 
ASCII file, convert it to binary and then place the message in 
a structure where the scaled and formatted data could be 
analyzed. Then GPSCRNCH.H "crunched" the numbers and computed 
the maximum, mean and standard deviation for the various time 


delays. 
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V. VERIFICATION OF OPERATIONAL CAPABILITIES 


GPS is a complicated new tool for navigators. Each 
receiver manufacturer employs unique and innovative techniques 
to improve its accuracy and flexibility. Motorola provides 
specifications and operational capabilities for its receiver 
and users have published various accuracies achieved using 
DGPS. It is prudent to verify those features important to this 
project’s success. Those two capabilities are examined below. 
They are: 


1. Verification of the manufacturer’s reacquisition time 
specifications. 


2. Determination of the effects of pseudorange correction 
delay in getting from the stationary receiver to the mobile 
receiver. 

In the following sections, these concerns are investigated. 


For each case, the experimental setup is described and the 


results are presented. 


A. REACQUISITION TIME VERIFICATION 

This project intends to use DGPS with one of the two 
antennae on a UAV. During some flight maneuvers, it is 
expected that the GPS antenna will be masked from one or more 
of the satellites it is using to calculate its position. It is 
important for the designers of the Kalman Filter to understand 


the normal reacguisition times given such a circumstance. In 


feo, 


Table IV in Chapter IV, Motorola presents its results for the 
reacquisition time of the PVT-6 receiver given the antenna was 
masked for 15, 30, 45 and 60 seconds. The method they used to 
collect this data is not well documented, and it is not 
indicated whether or not the receiver was in three dimensional 
fix mode or two dimensional fix mode. The difference could be 
Significant because the former requires four satellites fora 
fix and the latter only three. 

The reacquisition time was investigated using this 
project’s PVT-6 receiver. Table V shows the results next to 
the Motorola claims. Data was collected using the software 
supplied by Motorola with an evaluation kit. The antenna was 
completely masked 360° around above and below using four to 
Six layers of tin foil. The mask was left in place for the 
specified amount of time and then removed for at least 30 
seconds. Data files were imported into a spread sheet where 
the time marked data was analyzed. 

The results do not exactly agree with those that Motorola 
publishes, but they are close. Significantly, data obtained 
during daylight hours were within one standard deviation of 
Motorola’s claimed reacquisition times whereas data collected 
after daylight hours were not. No effort was made to confirm 
this suggested phenomenon nor was any effort made to correlate 
it to ionospheric thickness or tropospheric conditions such as 
rain or cloudiness. However, until the IMU is integrated with 


GPS using a Kalman Filter and the effects of a two to three 
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second delay in reacquisition time can be determined, 
investigation of this discovery is left for further study. 


Table V COMPARISON OF MOTOROLA AND EMPIRICAL REACQUISITION 
TIMES. 


a = 


Time Obscured Motorola l Standard 
Reacquisition Reacquisition Deviation 
Time (sec) Time (sec) (sec) 


13.754 





nO ..6 )2* 


de 7 he 


gree i 





* These observations were taken from 1230 to 1500. 
*x** These observations were taken from 2200 to 2400. 


B. DGPS ACCURACY AND PSUEDORANGE CORRECTION LATENCY 


The accuracy that DGPS can achieve is well documented. The 
PVT-6 is designed and built to permit DGPS. In fact, the 
Motorola Binary format includes master/slave DGPS commands and 
automatically incorporates the data in its calculations. It is 
worthwhile to validate the expected results especially in 
light of the fact the communications software to be used is 
not written by the manufacturers and may adversely affect the 
aceuracy . 

In designing the test for DGPS, the length of the baseline 
between the two receivers waS an issue that was researched. 
The results indicated the baseline range between the two 


receivers did not significantly affect DGPS accuracy up to 


oi 


about 50 km. Beyond 50 km, out to about 1000 km, the effects 
slowly increased the error (Kremer, 1990, pp. 307-312). This 
project’s maximum baseline is not expected to exceed 15 km. 
For all intents and purposes, a zero meter base line is as 
good as a 1000 meter baseline and certainly is easier to 
test. Therefore, this test simply mounted the two antennae 
side by side at a pre-surveyed location (See Appendix F). 
Data was collected by sending pseudorange corrections from the 
master antenna to the slave antenna using the software driver 
in Appendix E. The PVT-6 puts out a position message at a rate 
of once per second. These positions were differenced using the 
ECEF equations in Appendix A, and the mean, maximum and 
standard deviation were calculated. 

Table VI shows the DGPS accuracy with various delays in 
transmitting the correction to the slave receiver and Figure 
39 graphically illustrates the same results. The delay was 
artificially introduced through the function call to 
"Pr delay" in the header file GPSDELAY.H. The delays were in 
integer amounts of seconds as shown below. 

Interestingly, Motorola technicians quote adverse effects 
of pseudorange corrections received after 90 seconds. The 
results above suggest an immediate degrading effect of any 
amount of delay and that the error grows with delay time. It 
is hypothesized that the decrease in error at the 15 second 
delay point is a spurious result caused by any number of 


factors. It is left for further study to confirm this. 
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Table VI OBSERVATION OF PSEUDORANGE CORRECTION INPUT DELAY 
Om THe EFFECTIVE ACCURACY OF DIFFERENTIAL GPS. 


Average | Maximum | Standard 


Offset 
(m) 
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Standard Deviation 





seconds 
Figure 39 Graphical Representation Of The Data In Table VI 
Showing The Effects Of Time On The Accuracy Of DGPS 
Corrections. 
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VI. CONCLUSION 


This thesis reviews the major contemporary electronic 
navigation systems available and demonstrates that GPS offers 
the best accuracies and world wide coverage on a realtime 
basis. The user segment of GPS is explained, and its 
capabilities and deficiencies discussed. Differential GPS is 
presented and shown to be able to achieve enhanced accuracies 
sufficient to significantly contribute to this project. 

An extensive review of available GPS receivers was 
conducted and Motorola’s PVT-6 was found to combine all the 
benefits of GPS in a small, lightweight, low power and 
inexpensive package. Motorola’s proprietary binary interface 
protocol allows the user great flexibility with little 
sacrifice of convenience. In addition, the standardized RS-232 
serial communications port enables easy access to a variety of 
different systems. All things considered, it is the best 
choice over all other GPS receivers for this project. 

The PVT-6’s reacquisition times compare favorably with the 
manufacturer’s claims. When two of these receivers were used 
in a DGPS mode with a zero baseline, the results easily 
confirmed nominal results of one to three meters desired 
accuracy. The most encouraging result is the relative immunity 
of DGPS accuracy to short delays (one to seven seconds) of the 


pseudorange correction inputs. 
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The success of this project’s ultimate goal: to perform an 


autonomous landing, will depend on improving processing speeds 


and 


GPS 


accuracies. A few areas where further investigation into 
could contribute to this end are listed below. 


Conduct dynamic tests to confirm expected nominal 
accuracies using DGPS. 


Conduct additional data collection to determine if the 
accuracy of DGPS or reacquisition time is affected by the 
time of day or by the weather. 


Collect more data to determine if psuedorange correction 
latency steadily degrades or has anomalies (wherein the 
distance error actually improves with time) as Figure 39 
Suggests at the 15 second mark. 


Determine the effects of the aircraft’s RF environment on 
the GPS performance. 


Investigate whether programming in assembly or some other 
language might improve processing speeds. 


Determine if the PVT-6’s capability to compute two 
dimensional position by constraining altitude (relying 
solely on barometric altitude sensor data for altitude 
inputs) might improve robustness of GPS_ satellite 
trackinge 


Write a serial port driver that is independent of the PCL- 
744 multiple serial port I/O card in order to provide 
flexibility in using the PVT-6 and possibly reduce the 
total hardware weight through the elimination of the PCL- 
744 itself. 


Write a mission initialization interface routine to setup 
up the PVT-6 receivers. 


Write a post-processing routine to interpret and present 
the recorded results in a desired format. 


Write a real time graphical interface to plot the 
coordinates on a map. Also include either a separate 
Window or color coded routine to be able to see current 
altitude. 
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GPS has revolutionized electronic navigation. The airline 
industry and other nations decry the U.S. DOD’s introduction 
of Selective Availability. They are actively lobbying the U.S. 
government to force the DOD to turn it off and leave it off. 
It is an uphill battle but not a hopeless one. In time, as 
civilian applications become common place in our society, the 
DOD may acquiesce and turn it off. If that happens, GPS 
accuracies achieved in this project will undoubtedly improve 
and may warrant reconsideration in the way they are used in 


the Kalman Filter that integrates data from all the sources. 
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APPENDIX A: COORDINATE TRANSFORMATIONS 


A. GEODETIC COORDINATE TRANSFORMATION TO ECEF COORDINATES 
Given the geodetic coordinates longitude (A), latitude (¢) 
and altitude (h), the Earth Centered Earth Fixed coordinates 


(X,Y,Z) are (Leick, 1990, p. 184-185): 


X= (N+h) cosgd@ cosa 
Y=(N+h) cosgd sina 
Z={[{N(1-e7) +A] sing 


where h is the height above the ellipsoid and 


pee = sax 


Aa =s)p 
a 


f = 





a = semi-major axis 


b = semi-minor axis 
a 
v1 - e* sing 


B. ECEF COORDINATE TRANSFORMATION TO GEODETIC COORDINATES 
Given the ECEF coordinates (X,Y,Z), the Geodetic 


coordinates (A,¢,h) are (Leick, 1990, p. 184-185): 
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APPENDIX B: EPHEMERIS ALGORITHM 


The following parameters are provided continuously for 
each satellite of the GPS constellation. Their exact placement 
and description within the ephemeris message is defined by the 
Interface Control Document, ICD-GPS-200 as published by 
Rockwell International Space Systems Division. Any serious 
manipulation of the ephemeris algorithm using GPS data 
requires careful consideration of this document. There are 
several important, yet subtle, facts about the time tagging 


and the corresponding validity of the Issue of Data messages. 


My Mean Anomaly at Reference Time 

An Mean Motion Difference from Computed Value 

2 Eccentricity 

CA Square Root of the Semi-Major Axis 

NX Longitude of Ascending Node of Orbit Plane at 
Weekly Epoch 

a, Inclination Angle at Reference Time 

QW) Argument of Perigee 


OMEGADOT Rate of Right Ascension 
IDOT Rate of Inclination Angle 


c Amplitude of the Cosine Harmonic Correction 
Term to the Argument of Latitude 


Amplitude of the Sine Harmonic Correction Term 
to the Argument of Latitude 
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iw Amplitude of the Cosine Harmonic Correction 
Term to the Orbit Radius 


On Amplitude of the Sine Harmonic Correction Term 
to the Orbit Radius 


Ce Amplitude of the Cosine Harmonic Correction 
Term to the Angle of Inclination 


= Amplitude of the Cosine Harmonic Correction 
Term to the Angle of Inclination 


re Reference Time Ephemeris 

IODE Issue of Data (Ephemeris) 

The following tables list the number of bits, the location 
within words three through ten (omitting the parity bits), the 
scale factor of the LSB, the ranye and the units for each of 
subframes two (Table VII) and three (Table VIII). Parameters 
indicated with an asterisk are to be two’s complimented with 
the sign bit in the MSB. Unless otherwise listed in the value 
range column, effective range is the maximum range attainable 
with the indicated bit allocation and scale factor. 

The following elements of the ephemeris algorithm are 
listed in order of expected solution. All the required 
information is given or may be found using the above from the 
GPS ephemeris message. 

WGS 84 value of the earth’s 


universal gravitational 
sec? parameter 


2 
Men oec005 « 1074 Meters: 


wo at 


Table VII BIT ALLOCATION AND SCALING FACTORS FOR WORDS THREE 
THROUGH TEN OF SUBFRAME TWO. (ROCKWELL, 1987, PP. 60-81) 


Parameter| No. of Location |LSB scale Value Units 
bits of bits factor range 
IODE 1-8 see ICD- 
— 200 


>i ats: | 


An 16% 25-40 oo semi- 
circles/ 

sec 

My 32* 41-72 ype semi- 
err. | 





eS _. | 73-88 88 | | radians” 


peeesi 
/ 421-136 136 pa eeeGliens 
137-168 —— 
it, | 16 | 169-184 pes 604,784 


185-192 | lsee ICD-| 
GPS-200 | 





*Parameters indicated are signed integers. All others are unsigned. 


WGS 84 value of the earth’s 





Q = 7.2921151467 x 10° rau rotation rate 
| (JA)? Semi-major axis 
Computed mean motion 
n) = eS rad 
A? sec 
caer Time from ephemeris reference 
kK oe epoch 


In the above equation, t is GPS system time at time of 
transmission (i.e., GPS time corrected for transit time 
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Table VIII BIT ALLOCATION AND SCALING FACTORS FOR WORDS 
THREE THROUGH TEN OF SUBFRAME THREE. (ROCKWELL, 1987, PP. 
60-81) 


Parameter Bit LSB scale Value 
location factor range 
SS radians 
32* 17-48 pe semi- 
circles 
| ase | ao-ca | 2” | | radians | 


Be 65-96 2" semi- 
rey 


es | 97-112 | -112 eee 


32* 113-144 pe semi- 
circles 
OMEGADOT 145-168 semi- 
circles/ 
Sec 
oe 169-176 see ICD- 
GPS-200 
IDOT 177-190 semi- 
circles/ 
sec 


other 194-192 see. LEDs 
GPS-200 





*Parameters indicated are signed integers. All others are unsigned. 


range/speed of light . Furthermore, t, shall be the actual 
total time difference between the time t and the epoch time 
t,, and must account for beginning or end of week crossovers. 
That is, if t, is greater than 302,400 seconds, subtract 
604,800 seconds from t,. If t, 1s less than 302,400 seconds, 
add 604,800 seconds to t,. 


Corrected mean motion 


Mean anomaly 


The next equation solves for the eccentric anomaly (E,) by 
iteration. M, and e are known from above. 


M, 
Let E> 
i - e 
Then M,=E,-esinE ff MM, *M, then E, = E, 
otherwlse repeat and 
M, - M, compare again. 
A =F ee ee A 
ji oc colons (Tattersfield, 1984, p. 36) 
E, = E, +4, 


Having found E,, the rest of the solution is simply solving 
equations. 


Ji - e? sinE, True anomaly 
Vv, = arctan) ————_———— 
COS, == 56 
a Argument of latitude 


Second harmonic perturbations: 


du, =C,, sin2®, + C,, cos2®, Argument of lat correction 
6x, =C,. cos2%, + C,, sin2®, Radius correction 
61,=C;..cos2®, + C,. sin2@, Correction to inclinateen 


Corrected argument of latitude 


Corrected radius 


i, _ i, ‘ Si, + (IDOT) t, Corrected inclination 
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- b COSsU, 


Yo = LC, Sinu, 


ee, / ; ; 

X= x COS hile — y;. COS1, Sint, 
/ : / 

Prue Sima, — Y,COSl, COSA) 


_ 7 . . 
me — Y, Sini, 


Positions in orbital plane 


Corrected longitude of 
ascending node 


Earth-Centered, Earth-Fixed 
(ECEF) coordinates 
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Antenna 


Antenna to Recetver 
Interconnection 


Serial Output 


Accuracy 


Altitude 


Operating Temperature 


Humidity 


Physical Dimensions 


Weight 


Switched Power 


"Keep-Alive® BATT Power 


Power Consumption (typical) 


Receiver Interconnection 


MTBF (Mean Time Before Failure) 


APPENDIX C: PVT-6 SPECIFICATIONS 


= 


Peay Big ve ibabty a5 bee. 4: 2k 
[Sl as ‘ 


— te 


wes SRE Se ie ¢ Pe a 
By % 5 us a abe ti ae ets sve oleae “ Descriptior oy ae a aaraa She i eee ae a - sey | 
a Py <A 
[Me 


arieta ie 


ancy 1000 Knots (515 m/sec) 
Acceleration 4g 
Jerk 5mvs° 


Active microstrip patch Antenna Module 
Powered by Receiver Module (5 Vdc @ 25 mA) 


Single coaxial cable (6dB max loss at L1; 1575.42 MHz) 
Typical length with RG58 coaxial cable; 20 feet (6 meters) 


RS232C Interface 


Less than 25 meters, SEP (without SA) 
* DoD may invoke Selective Availability (SA), potentially Coe 
accuracy to 100 meters (2d RMS) 


60,000 feet (18 Km) (maximum) 


Active Antenna ~40°C to +100°C (-40°F to +212°F) 
Receiver Module -30°C to +80°C (-22°F to 176°F) 


95% non-condensing +30°C to +60°C (+86°F to +140°F) 


Receiver PCB 3.94 X 2.76 X 0.65 in (100 X 70 X 16.5 mm) 
Antenna Module 4.01 (dia) X 0.89 in (102 (dia) X 22.6 mm) 


Receiver and Housing 4.5 oz (128 g) 
Antenna Module 4.8 02 (136.2 g) 


9 to 16 Vdc or 
5 +0.25 Vde 


| 


4.75 - 16 Vde; 0.3 mA 


1.3 W @ 5 Vde input 
1.8 W @ 12 Vdc input 


AMP 10 pin, #104369-4 on Receiver Module 
- RS232C 
- Power input 

OSX connector on Receiver and Antenna 


55,000 hours (estimated) 
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APPENDIX E: SOFTWARE DRIVER SOURCE 


A. GPSAN.C 


GPSAN.C 


Reads ascii position data "@@Ba" from a file, converts it 
it to binary and places the scaled and formatted messages 
in an array of structures. 


Input is the name of the ascii file without the extension. 
Output is a file with the analyzed data from GPSCRNCH.H. 


#include <string.h> 

#include <stdlib.h> 

#include <stdio.h> 

#include "“c:\borlandc\twite\gpstruct.h” 
#include "c:\borlandc\twite\gpsconv.h" 
#include "c:\borlandc\twite\gpscrnch.h" 
#include "c:\borlandc\twite\gpsfun.h" 


void main(void) 


#define COL 130 
#define NEWCOL 66 


char name[40], path[30] = "c:\\borlandc\\etwork\\"; 

int i,j,kK,m,n,s,master,slave, ROW; 

ONEBYTE first, second, third, inbuff[COL]; 

ONEBYTE **mbuff, **sbuff, **mbuffbin, * *sbuffbin; 
struct T POS CHAN STATUS “*marst, *sarst, pos_chan; 
FILE *fpname; 


printf("\nEnter the input filename (without the extension(??.ext): "); 
scanf(" %s",name); 

strcat(path, name); 

strcat(path, ".txt"); 

printf("%s",path); 

fpname = fopen(path,"r"); 


printf("\nEnter the number of rows (individual messages): "); 


scanf(" %d", &ROW); 
printf("\n"); 
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CODE 


/* Create mbuff and sbuff. 

oe 

mbuff = calloc( ROW, sizeof{( ONEBYTE * )); 
sbuff = calloc( ROW, sizeof( ONEBYTE * }); 
for (i=0; i< ROW; i+ +) 


mbuffli] = calloc( COL, sizeof{( ONEBYTE )); 
sbuffli} = calloc( COL, sizeof{ ONEBYTE )); 


}; 
while (1) 
{ 
/* Read in file to ascii arrays mbuff and sbuff. 
ay, 
if(k >= 70) 
break; 
first = fgetc(fpname); 
if ( first = = '@' ) 
{ 
second = fgetc(fpname); 
if (second = = ‘@’ ) 
{ 
third = fgetc(fpname); 
/* Store inbuff array in either mbuff or sbuff. 
af 
if (third == '@") 
fread( sbuffl[s], sizeof(ONEBYTE), COL, fpname ); 
kK = 0; 
s+ +: 
} 
else 
mbuff(mJ[O} = third; 
for (i=1; i1<COL; i+ +) 
mbuff(m][i} = fgetc(fpname); 
1 8 
m+ +; 
}: 
a; 
}; 
k++; 


}; /* end while */ 
fclose(fpname); 


master = m; 
slave = sS; 


<)e) 


/* Create mbuffbin and sbuffbin. 
Td 
mbuffbin = calloc( master, sizeof( ONEBYTE * )); 
for (i=0O; i< master; i+ +) 
mbuffbin[i] = calloc( NEWCOL, sizeof( ONEBYTE )); 
sbuffbin = calloc( slave, sizeof( ONEBYTE ®* )); 
for (i=0; i<slave; 1+ +) 
sbuffbin[i] = calloc( NEWCOL, sizeof( ONEBYTE )); 


/* Convert hex buffer to binary buffer. 
i: 
for (i = O; i < master; i+ + ) 


mbuffbin[iJ[O] = mbuff[iJ[O]; 
mbuffbin{i][1] mbuff[iJ[1]; 
k = 2; 

fort) = 27j'— COL 4 —Z) 


mbuffbin[iJ[k] = (Hex(mbuffli]fj]) << 4) | Hex(mbuffli}[{j + 1)); 
kK += 1; 
}; 
}; 


for (i = O: i < slave; i+ + ) 


sbuffbin[iJ[O] = sbuffli][O}; 
sbuffbin[i}[1] = sbuff[il[1]; 
Case 

for {.j =22; = 26 Ok j- —Z) 


{ 

sbuffbin[iJ[k] = (Hex(sbuff[i}[j]) << 4) | Hex(sbufflifj+ 1)); 
Ko = 1 

hi 
} 


/* Free mbuff and sbuff memory. 
ey 
for (1\=0; i< ROW; i+ +) 


free(mbuffl[i]); 
free(sbuff[1]); 
\ 
free(mbuff); 
free(sbuff); 


/* Temporarily store mbuffbin in a file "Mtemp.bin" and free memory. 
7) 
fpname = fopen( "mtemp.bin", "wb" ); 
for (i\=O; i< master; i+ +) 
fwrite(mbuffbinl[i], sizeof(ONEBYTE), NEWCOL, fpname); 
for (1\=O; i< master; 1+ +) 
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free(mbuffbinl[i]); 
free(mbuffbin); 
fclose(fpname); 


/* Temporarily store sbuffbin in a file "stemp.bin" and free memory. 
if 
fpname = fopen( "stemp.bin", "wb" ); 
for (i'=O; i<slave; i+ +) 
fwrite(sbuffbin[i], sizeof(ONEBYTE), NEWCOL, fpname); 
for (i =O; i<slave; i+ +) 
free(sbuffbinl[i]); 
free(sbuffbin); 
fclose(fpname); 


/* Create an array of structures. 

ot | 

marst = calloc( master, sizeof( struct T POS CHAN STATUS )); 
sarst = calloc( slave , sizeof( struct T POS CHAN STATUS )); 


fpname = fopen( "mtemp.bin", "rb"); 
for (i = O; 1 < master; i+ +) 


fread(inbuff, sizeof(ONEBYTE), NEWCOL, fpname); 
Convert_Pos Chan_Status( inbuff ); 
marstli] = pos_chan; 

n 

fclose( foname ); 

remove("mtemp.bin"); 

fpname = fopen( "stemp.bin", "rb"); 

for (i = 0; 1 < slave; i+ +) 


fread(inbuff, sizeof(ONEBYTE), NEWCOL, fpname); 
Convert_Pos_Chan_Status( inbuff ); 
sarstli] = pos chan; 


}; 
fclose( fpname ); 
remove("stemp.bin"); 


Ave _Max_Stdev(marst, sarst, master); 


}; /* end main */ 
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B. GPSCONV.H 


/* GPSCONV.H 


Convert_Ephemeris() 

Convert_Pos_Chan_Status() 

Convert_xDOP() 

Convert_Sat_Range_Data() 

Convert Pos Chan_Status_Ext() 
ma 


#ifndef _GPSCONV_H 
#define _GPSCONV_H 


#include <math.h> 


#ifndef _GPSTRUCT H 
#include "“c:\borlandc\twite\gpstruct.h” 
#endif 


Lee ee eee eee ee ee a en ee Oe ee ee 


Convert_Ephemeris() 


This procedure converts the raw binary ephemeris data 
into a representation usable for post-processing 
applications. The equations coded in this procedure 

are taken from the ICD-GPS-200 (the GPS NAVSTAR ICD). 


Input ts an array with the binary ephemeris message. 
Although it is unscaled, it is in the array 
format of the [CD-GPS-200 format. 
Output is a structure with the scaled and formatted 
ephemeris message. 
ae 


void Convert_Ephemeris( ONEBYTE hex_ephemeris{]{3] ) 
{ 
/* Sign extends a neg hex byte into its TWOBYTE (integer) equivalent 
ey 
#define MACRO_SET_UPPER(word_to set) \ 
if ( word to set & 0x0080) word _to set ; = OxFFOO ; 


/* This MACRO takes two hex bytes and merges them into a single 
TWOBYTE (integer) value. 
ba 
#define MACRO MAKE EPH WORD(e_ wd, wil, b1, w2, b2) \ 
e wd = ((( UNSIGNED_TWOBYTE ) hex_ephemeris[w1][b1] << 8) | \ 
( UNSIGNED _TWOBYTE ) hex_ephemeris[w2][b2] ); 
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UNSIGNED _TWOBYTE word1, word2; /* Temporary storage for words built. */ 
struct T FLOAT_EPHEMERIS “*eph ; 


eph = malloc(sizeof( struct T FLOAT EPHEMERIS) ); 


/* Time of ephemeris 

a 

MACRO_MAKE_EPH_WORD(word1, 15, 0, 15, 1 ); 

eph->toe = (double) ( ((UNSIGNED FOURBYTE) word1) * 16); 





/* Mean anomaly at reference time ( m0 ). 
| */ 
| 
MACRO _MAKE_EPH_WORD(word1, 9, 2, 10, 0); 
MACRO _MAKE_EPH_ WORD(word2, 10, 1, 10, 2 ); 
eph->m0O = ( ( ( FOURBYTE ) word] << 16) | word2 ) * PITWOM31; 


/* Mean motion difference from computed value ( delta_n ). 
ba 

MACRO MAKE _EPH_WORD(word1, 9, 0, 9, 1 ); 
eph->delta n = ( TWOBYTE ) word] * PITWOM43; 


—-__-___—-—.- 


we 

MACRO _MAKE_EPH WORD(word1, 11, 2, 12, 0); 
MACRO_MAKE_EPH_WORD(word2, 12, 1, 12, 2 ); 

eph->e = ((( UNSIGNED FOURBYTE ) word1 << 16) | word2 }) * TWOM33; 


/* Eccentricity (e). 


/* Square root of the semi-major axis ( sqrt_a ). 

“y 

MACRO_MAKE_ EPH WORD(word1, 13, 2, 14, 0); 

MACRO. MAKE EPH WORD(word2, 14, 1, 14, 2); 

eph->sqrt_a = ( ( ( UNSIGNED FOURBYTE ) word1 << 16) | word2 ) * TWOM19; 


/* Longitude of ascending mode of orbit plane at 
weekly epoch ( omega 0 ). 
ry, 
MACRO_MAKE EPH WORD(word1, 16, 2, 17, 0); 
MACRO MAKE EPH WORD(word2, 17, 1, 17, 2); 
eph->omega 0 = (( ( FOURBYTE ) word] << 16) | word2 ) * PITWOM31; 


/* Inclination angle at reference time ( 10 ). 

oy 

MACRO _MAKE_EPH WORD(word1, 18, 2, 19, 0); 

MACRO MAKE EPH WORD(word2, 19, 1, 19, 2); 

eph->i0 = ( ( ( FOURBYTE ) word1 << 16) | word2 ) * PITWOM31; 


/* Argument of perigee ( w ). 

a 

MACRO MAKE EPH WORD(word1, 20, 2, 21, 0); 

MACRO MAKE EPH WORD(word2, 21, 1, 21, 2 ); 

eph->w = ((( FOURBYTE ) word1 << 16)! word2 ) * PITWOM31; 
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/* Rate of right ascension ( omega_dot ). 

| 

word] = ( UNSIGNED TWOBYTE ) hex ephemeris[22}[0}; 
MACRO SET UPPER(word1); /* Set upper bits if negative. */ 

MACRO MAKE _EPH WORD(word2, 22, 1, 22, 2); 

eph->omega dot = ((( FOURBYTE ) word] << 16) | word2 ) * PITWOM43; 


/* Rate of inclination angle (i dot ). 

ay 

MACRO MAKE EPH WORD(word1, 23, 1, 23, 2); 

word] = word! >> 2; 

if ( word] & 0x2000 ) /* Set upper bits if negative. */ 
word] | = OxCOO0O; 

eph->i dot = ( TWOBYTE ) word1 * PITWOM43; 


/* Amplitude of the cosine harmonic correction term to the 
argument of latitude (cuc). 

Gil 

MACRO MAKE EPH WORD(word1, 11, 0, 11, 1 ); 

eph->cuc = ( TWOBYTE ) word1 * TWOM29; 


/* Amplitude of the sine harmonic correction term to the 
argument of latitude (cus). 

ad 

MACRO_MAKE EPH WORD(word1, 13, 0, 13, 1 ); 

eph->cus = ( TWOBYTE ) word1 * TWOM29; 


/* Amplitude of the cosine harmonic correction term to the 
orbit radius (crc). 

| 

MACRO MAKE EPH WORD(word1, 20, 0, 20, 1 ); 

eph->crc = ( TWOBYTE ) word1 * TWOMS; 


/* Amplitude of the sine harmonic correction term to the 
orbit radius (crs). 

af 

MACRO MAKE EPH WORD(word1, 8, 1, 8, 2 ); 

eph->crs = ( TWOBYTE ) word] * TWOMS5; 


/* Amplitude of the cosine harmonic correction term to the 
angle of inclination (cic). 

ih 

MACRO MAKE EPH _WORD(word1, 16, O, 16, 1); 

eph->cic = ( TWOBYTE ) word! * TWOM29; 


/* Amplitude of the sine harmonic correction term to the 
angle of inclination (cis). 

bar 

MACRO MAKE EPH WORD(word1, 18, 0, 18, 1); 

eph->cis = ( TWOBYTE ) word] * TWOM29; 
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/* Estimated group delay differential ( tgd ). 

= 

word1 = ( UNSIGNED_TWOBYTE ) hex_ephemeris(4][{2]; 
MACRO_SET_UPPER(word1); /* Set upper bits if negative. */ 
eph->tgd = ( TWOBYTE ) word1 * TWOM31; 


/* Clock data reference time ( toc ). 

oa 

MACRO MAKE_EPH_WORD(word1, 5, 1, 5, 2); 
eph->toc = word! * 16.0; 


/* Polynomial coefficient ( af2 ). 

a. 

word] = ( UNSIGNED_TWOBYTE ) hex_ephemeris(6)[0]; 
MACRO SET_UPPER(word1); /* Set upper bits if negative. */ 
eph->af2 = ( TWOBYTE ) word! * TWOM55; 


/* Polynomial coefficient ( af1 }. 

a 

MACRO MAKE EPH WORD(word1, 6, 1, 6, 2); 
eph->afl = ( TWOBYTE ) word! * TWOM43; 


/* Polynomial coefficient ( afO ). 

| 

word] = ( UNSIGNED _TWOBYTE ) hex _ephemeris([7}(0}; 
MACRO_SET_UPPER(word1); /* Set upper bits if negative. */ 

MACRO MAKE EPH WORD(word2, 7, 1, 7, 2); 

word2 = word2 & OxFFFC ; /* clear 2 bits */ 

eph->afO = ((( FOURBYTE }) word1 << 16) | word2 ) * TWOM3125; 


/* Semi-major axis (a }. 
ay 
eph->a = eph->sqrt_a * eph->sart_a; 


/* Corrected mean motion (n ). 
“) 
eph->n = SQRMU / ( eph->a * eph->sqrt_a) + eph->delta_n; 


/* Square root of eccentricity*2 subtracted from 1 ( esart ). 
Be dart = sqrt (1.0- eph->e * eph->e ); 

/* Omega dot - WE ( wk1 ). 

eee = eph->omega dot - WE; 

/* OmegaO - WE * toe ( wkO }). 

etc = eph->omega 0 - WE * eph->toe; 


/* Sine w ( sin_w ). 
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a] 
eph->sin_w = sin ( eph->w ),; 


/* Cosine w (cos w). 
wy) 
eph->cos_w = cos ( eph->w ); 


/* Square root of semi-major axis * eccentricity * FCONST. 
| 
eph->sqrta_e fconst = eph->sqrt_a * eph->e * FCONST; 


/* Issue of data, ephemeris. 
“ 
eph->iode = hex ephemeris[23][0]; 


return *eph; 


}; 


ee ee ee ee Ee ee om ie eee in gece ee ore 


Convert _Pos_Chan_Status() 


This procedure parses and decodes a position/status/data 
record received from the Motorola PVT-6. The flow of this 
procedure is as follows : 

- pass in the BINARY buffer containing 66 bytes 

- read the BINARY data from the buffer and convert directly 

into the final data structure. 


Input is a buffer of type ONEBYTE with a single binary 
position/status/data message. 
Output is the structure with the scaled and formatted 
position/status/data message. 
7 


struct T POS CHAN STATUS Convert_Pos_Chan_Status(ONEBYTE posl]) 


{ 

struct T POS CHAN_STATUS  pos_chan; 
UNSIGNED ONEBYTE i, J; 

UNSIGNED ONEBYTE tempchar ; 
UNSIGNED FOURBYTE tempu4byte ; 
FOURBYTE temps4byte ; 

UNSIGNED FOURBYTE nsecs ; 

double degrees, minutes, seconds ; 


/* Read and scale the of the data. 
a 

pos chan.month = pos[(2] ; 

pos chan.day = pos([3] ; 


tempchar = pos[4] ; 
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pos chan.year = ( tempchar << 8) + pos([5] ; 


pos_chan.hours = pos[6] ; 
pos chan.minutes = pos[7] ; 


tempchar = pos[8]; /* integer seconds */ 
tempu4byte = pos[9] ; 

tempu4byte = (tempu4byte << 8) + pos[10]; 
tempu4byte = (tempu4byte << 8) + pos[{11]; 
tempu4byte = ( tempu4byte << 8) + pos[12] ; 


pos chan.seconds = (double) tempchar + \ 
( ( (double) tempu4byte ) / 1.0E +9 ); 


temps4byte = pos[13] ; 

temps4byte = (temps4byte << 8) + pos[14]; 
temps4byte = ( temps4byte << 8) + pos[(15]; 
temps4byte = ( temps4byte << 8) + pos[16] ; 
degrees = (double) temps4byte * MSECS TO DEGREES ; 


II 


pos chan.latitude.degrees = (TWOBYTE) degrees ; 

if (degrees < 0) 

degrees = fabs ( degrees ) ; 

minutes = (degrees - (TWOBYTE) degrees ) * 60.0; 

pos _chan.latitude.minutes = (TWOBYTE) ( minutes ) ; 

pos _chan.latitude.seconds = ( minutes - (TWOBYTE) minutes ) * 60.0 ; 


temps4byte = pos[17]; 

temps4byte = (temps4byte << 8) + pos[18]; 
temps4byte = ( temps4byte << 8) + pos[19]; 
temps4byte = ( temps4byte << 8) + pos[20] ; 
degrees = (double) temps4byte * MSECS TO DEGREES ; 


pos chan.longitude.degrees = (TWOBYTE) degrees ; 

if ( degrees < O) 

degrees = fabs ( degrees ) ; 

minutes = ( degrees - (TWOBYTE) degrees ) * 60.0 ; 

pos chan.longitude.minutes = (TWOBYTE) ( minutes ) ; 

pos _chan.longitude.seconds = ( minutes - (TWOBYTE) minutes ) * 60.0 ; 


temps4byte = pos[21]; 

temps4byte = ( temps4byte << 8) + pos[22] ; 
temps4byte = ( temps4byte << 8) + pos[23] ; 
temps4byte = (temps4byte << 8) + pos[24] ; 
pos_chan.datum_ height = (double) temps4byte / 100.0 ; 


temps4byte = pos[25] ; 


temps4byte = ( temps4byte << 8) + pos[26] ; 
temps4byte = ( temps4byte << 8) + pos(27] ; 
temps4byte = ( temps4byte << 8) + pos[28] ; 


pos chan.msl_height = (double) temps4byte / 100.0 ; 


LO 


tempchar = pos[29] ; 
pos chan.velocity = (double) ( ( tempchar << 8) + pos[30] ) / 100.0 ; 


tempchar = pos([31] ; 
pos chan.heading = (double) ( ( tempchar << 8) + pos[32])/ 10.0; 


tempchar = pos(33] ; 
pos_chan.current_dop 


ll 


(double) ( ( tempchar << 8) + pos[34] ) / 10.0; 


pos chan.dop_ type = pos[35] ; 
pos chan.visible sats = pos[36] ; 
pos chan.sats tracked = pos([37] ; 


j=0; 

for (i = 0; i < NUM CHANNELS; i+ +) 
{ 
pos _chan.channelli].svid = pos[38 + j); 
pos chan.channelli].mode = pos[39 + j); 
pos_chan.channelli].strength = pos[40 + j]; 
pos_chan.channelli].flags = pos[41 + j); 
j+ =4:; 
} 


pos chan.rcvr_ status = pos[62]; 


return pos chan; 


jie 


#endif /* GPSCONV_H */ 
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ae GPSCRNCH.H 


eeorockNCH.A 


Used to crunch the position data to determine the mean, 
maximum, and standard deviation of the position error 
obtained by DGPS. Geodetic coordinates are converted to 
ECEF coordinates and the distance formula is used to compute 
the distance error. 

Input is an array of structures that contain the scaled 

and formatted position/chan/status message data. 
The user is asked to supply the output file name 
without the extension. 

Output is a file with the distance error for each data 
point, and the mean, max and standard deviation of 
all the data points. 

| 


#ifndef _GPSCRNCH H 


#define _GPSCRNCH _H 


#include <math.h> 
#include <stdio.h> 
#include <string.h> 


#ifndef GPSTRUCT H 
#include “c:\borlandc\twite\gpstruct.h” 
#endif 


FILE “{p: 


void Ave Max Stdev(struct T POS CHAN STATUS “mast, \ 
struct T POS CHAN STATUS “*slav, int master) 
{ 


int i, count =O, start; 

char name[40], outpath[30] ="c:\\borlandc\\twite\\"; 

double mlat, mlon, malt, slat, slon, salt, MN, sN, mean, distsum =O; 
double mx, my, mz, sx, sy, sz, dist{200], max =-1.0, stdev, distsum2 =O; 


printf("Enter the output filename without the extension (???.txt): "); 
scanf("%s",name); 

strcat(outpath, name); 

strcat(outpath, ".out"); 

fp = fopen(outpath,"w"); 

fprintf(fp, "\n%s\n" ,outpath); 


printf("Enter the starting data point number: "); 
scanf("%d",&start); 
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for (i\=Sstart; i< master; i+ +) 
{ 
/* Test to ensure the seconds match */ 
if ( (ONEBYTE) mastli].seconds = = (ONEBYTE) slav[i].seconds ) 


/* Convert to mx, my, mz, Sx, SY, SZ. 

i 

miat = (mast[i].latitude.seconds * (1/60) + \ 
mast[i].latitude.minutes) * (1/60) + mastli].latitude.degrees; 

mlon = (mast[i].longitude.seconds * (1/60) + \ 
mast[i].longitude.minutes) * (1/60) +mast[i].longitude.degrees; 

malt = mast(t].datum_height; 


slat = (slav[i].latitude.seconds * (1/60) + \ 
Slav[i].latitude.minutes) * (1/60) +slavli].latitude.degrees; 
slon = (slav[i].longitude.seconds * (1/60) + \ 
slav[i].longitude.minutes) * (1/60) +slavli].longitude.degrees; 
salt = slavli].datum_height; 


mN = SEMIMAJOR_ AXIS / sqrt( 1 - E SQUARED * sin(mlat) ); 
sN = SEMIMAJOR AXIS / sqrt( 1 - E SQUARED * sin(slat) ); 


mx = (mN + malt) * cos(mlat) * cos(mlon); 
my = (mN + malt) * cos(mlat) * sin(mlon); 
mz = (mN * (1-E SQUARED ) + malt) * sin(miat); 


sx = (SN + salt) * cos(slat) * cos(slon); 
sy = (sN + salt) * cos(slat) * sin(slon); 
sz = (sN *(1-E SQUARED ) + salt) * sin(slat); 


distli] = sqrt( (mx-sx)*(mx-sx) + (my-sy)*(my-sy) + (mz-sz)*(mz-sz) ); 
if ( distli] > 1000 ) 


forintf(fp,"\nMtime = %2d:%10.7f <-> %2d:%10.7f = Stime  \ 
Dist #%3d = XXXXX",mastli].minutes, mast[i].seconds, \ 
slav[i].minutes, slav[i].seconds, i); 

continue; 


iF 


/* Compute average, max and standard deviation. 
i 

count+ +; 

if ( distli] > max ) max = distil]; 

distsum + = distlil; 

distsum2 + = dist{i] * dist[i]; 


/* Print out the results to a file. 

a 

forintf(fp, "\nMtime = %2d:%10.7f <-> %2d:%10./7f = Stime  \ 
Dist #%3d = %5.2f",mastli].minutes, mastli].seconds, \ 


Tne 


slavli].minutes, slavli].seconds, i, distli]); 


} 


else 


{ 

fprintf(fp,"\nMtime = %2d:%10.7f <-> %2d:%10.7f = Stime \ 
Dist #%3d = XXXXX",mastli].minutes, mastli].seconds, \ 
slavli].minutes, slav[i].seconds, i); 


/* Compute the mean and standard deviation. 

nt, 

mean = distsum / count; 

stdev = sart( (distsum2/count) - (mean * mean)); 


fprintfifp, "\n\nMean = %5.2f\nMax = %5.2f\nStdev = %5.2f\n",\ 
mean, max, stdev); 


fclose(fp); 
}; /* end Ave Max Stdev */ 


#endif /* GPSCRNCH H */ 


ga 


D. GPSDEFIN.H 


/* GPSDEFIN.H 


Store most of the definitions here. 
ie 
/ 


#ifndef GPSDEFIN H 
#define _GPSDEFIN_H 


[oe ae type definitions eo eager ey, 


typedef char ONEBYTE; /* data type allowing for the use of 8 bits. */ 
typedef short TWOBYTE; /* data type allowing for the use of 16 bits. */ 
typedef long FOURBYTE; /* data type allowing for the use of 32 bits. */ 


typedef unsigned char UNSIGNED ONEBYTE; 
typedef unsigned short UNSIGNED_TWOBYTE; 
typedef unsigned long UNSIGNED _FOURBYTE; 


/* unsigned 8 bits */ 
/* unsigned 16 bits */ 
/* unsigned 32 bits */ 


[prreeee*e** constant definitions  ~ / 


#define Pl 
#define WE 


3.1415926535898 
7.2921151467E-5 


/* GPS value of Pi ns 
/* WGS84 earth's rot rate (rads/sec) */ 


#define SQRMU  1.99649818432174E7 /* Square root of MU. is | 


#define FCONST -4.442807633E-10 /* FCONST is used in relativistic 

clock correction. Defined on 

page 73 of ICD-GPS-200. wy, 
#define INVERSE FLATTENING 298.257223563 /* WGS-84 unitless */ 
#define FLATTENING 1/INVERSE FLATTENING 
#define SEMIMAJOR_AXIS 6378137.0 /* WGS-84 meters */ 
#define E SQUARED 2*FLATTENING-FLATTENING *FLATTENING 


#define TWOM5 
#define TWOM11 
#define TWOM19 


0.03125000000000 
(TWOMS / 64.0) 
(TWOM11 / 256.0) 


(oS 


#define TWOM20 
#define TWOM21 
#define TWOM23 
#define TWOM24 
#define TWOM27 
#define TWOM29 
#define TWOM30 
#define TWOM31 
#define TWOM33 
#define TWOM38 
#define TWOM43 


(TWOM19 / 2.0) 
(TWOM20/ 2.0) 
(TWOM21 / 4.0) 
(TWOM23 / 2.0) 
(TWOM24 / 8.0) 
(TWOM23 / 64.0) 
(TWOM29 / 2.0) 
(TWOM30 / 2.0) 
(TWOM31 / 4.0) 
(TWOM33 / 32.0) 
(TWOM38 / 32.0) 


AeakZ. 


#define TWOM50 (TWOM43 / 128.0) 
#define TWOM55 (TWOM50 / 32.0) 
#define TWOM3125 (0.25 * TWOM31) 


#define TWOP11 2048.0 ie ate 7 
#define TWOP12 (TWOP11 * 2.0) 
#define TWOP14 (TWOP12 * 4.0) 
#define TWOP16 (TWOP14 * 4.0) 





#define TWOP12 INTEGER 4096 Bele ia Faia 
#define TWOP16_ INTEGER 65536 fad WOwwet 


| #define PITWOM31 (PI * TWOM31) 
#define PITWOM43 (PI! * TWOM43) 


aa Twite definitions el ala 


#define FALSE 0 

#define TRUE 1 

#define RX 0 /* receive */ 

#define TX 1 /* send */ 

#define MGPS_PORT 3 /* master gps receiver port */ 
#define Sore PORT 4 /* slave gps receiver port */ 


#define NUM CHANNELS 6 
#define MSECS TO DEGREES (( 1.0/ 1000.0 ) / 3600.0 ) 


#define EPH NUM WORDS 24 /* at 3 bytes/word */ 


#endif /* GPSDEFIN_H */ 
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E. 


he 


7. 


GPSDELAY.H 


GPSDELAY.H 


Used to introduce an artificial delay in transmitting 
psuedorange corrections from the master to the slave. 


PR_delay() accepts a pointer of type PANDL to the buffer 
containing the psuedorange correction message 
(@@Ce...) and an integer number of seconds it 
is to be delayed. The output is to the PCL-744 
serial 1/O card. 


#ifndef GPSDELAY_H 
#define GPSDELAY_H 


#include <stdio.h> 

#include <string.h> 

#include <stdlib.h> 

#include <time.h> 

#include <dos.h> 

#include “c:\borlandc\twite\head-c.h" 


#ifndef ©=GPSTRUCT _H 
#include "c:\borlandc\twite\gpstruct.h” 
#endif 


#define MAXDELAY 300 


static PANDL — stor[{MAXDELAY]; /* Ring buffer mf 
Static int num_waiting = 0; /* Number of messages in queue “*/ 
Static int next_stor = 0; /* Next storage location to be used */ 
static int next_xmit = 0; /* Next message to be transmitted */ 


void PR_delay( PANDL *Ce msg, int delay ) 


PANDL *ptr; 


/* Copy the contents of the buffer with the Ce msg to a 
new buffer for holding during the delay. 

= 

ptr = malloc( sizeof(PANDL) ); 

ptr->len = Ce _msg->len; 

ptr->ptr = calloc( ptr->len, sizeof(ONEBYTE) ); 

memcpy( ptr->ptr, Ce_msg->ptr, Ce _msg->len ); 


/* Store the structure of type PANDL in the ring buffer "stor[]". 


a 
if (num_waiting < MAXDELAY ) 
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{ 


Stor{next_stor] = “ptr; 
num_waiting+ + ; 

next stor+ + ; 

if (next_stor = = MAXDELAY) 


next_stor -= MAXDELAY; 
} 


else 


{ 

printf("There are %d \"@@Ce...\" messages in the ring buffer.\n",\ 
MAXDELAY); 

printf("One correction message was lost."); 

free(ptr- > ptr); 

free(ptr); 

return; 


} 


/* Transmit the oldest message if there are more messages 
than the delay length. The delay is premised upon the 
fact that correction messages come approximately one second 
apart. 

st | 

while ( (num_waiting - delay) > O ) 


{ 


sio_putb( SGPS_ PORT, stor[next_xmit].ptr, stor{next_xmit].len ); 


free(stor[next_xmit].ptr); 


next_xmit+ + ; 
if (next_xmit = = MAXDELAY) 


next_xmit -= MAXDELAY; 
num_waiting-- ; 


i 


free(ptr); 
return; 


} 
#Hendif /* GPSDELAY_H */ 


ED 


F. GPSFUN.H 


/* GPSEWN EH 
A collection of functions to manipulate GPS messages. 


Hex() 
Check_sum() 
Bin to_ascii() 
Msg _verification() 
Parse_inbuff() 
Master gps() 
Slave_gps() 

al 


#ifndef —GPSFUN_H 
#define —GPSFUN_H 


#include <stdlib.h> 
#include <string.h> 
#include <stdio.h> 
#include "c:\borlandc\twite\head-c.h” 


#ifndef ©GPSTRUCT _H 
#include "c:\borlandc\twite\gpstruct.h" 


#endif 
FILE *m_all_file; /* master file for all data */ 
|| BE *m_pos file; /* master file for pos data */ 


/* Global matrix for converting ephemeris 

ry 

UNSIGNED _ONEBYTE hex_ephemeris[24][3] ; 

/* [word][byte] : 24 words (3 bytes per word) corresponding to 
words 3-10 of subframes 1-3 of the satellite nav message. 

a 


ee eS Se te ee ee ee ny eee ee eee ee eee 


This procedure converts an ASCII value to its hex 
equivalent. 


Input is a single ONEBYTE. 


Output is a single ONEBYTE. 
WE 


1.6 








ONEBYTE Hex(ONEBYTE ASClivalue) 


{ 

/* O is Ox30(48) -> 9 is 0x39(57) 
A is 0x41(65) -> Fis 0x46(70) 
ais 0Ox61(97)->  f is Ox66(102) 

dl 


if ((ASClIlvalue) < 0x40) 
return((ASCllvalue) - 0x30); 
else 
if ((ASClIIlvalue) < Ox60) 
return((ASCIIvalue) - 0x37); 
else 
return((ASCllvalue) - 0x57); 
}; 
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Check_sum() 


Calculates the checksum of a message in a buffer of type 
ONEBYTE. Since the message is assumed to be in the Motorola 
Binary Proprietary format, all but the last three bytes 

are used in computing the checksum. 


Input is the pointer of type PANDL. 
Output is the checksum of type ONEBYTE. 
*/ 
ONEBYTE Check _sum( PANDL *msg ) 


{ 
ONEBYTE chk sum = O; 
int |; 


/* Compute the Exclusive-Or of all but the last three bytes 


ie 

for (1 = 0; 1 < (msg->len - 3); 1+ + ) 
chk_sum *= msg->ptrlil; 

return chk_sum; 

}; 
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Bin_to_ascii() 


Converts a message from binary to its hexidecimal 
ASCII equivalent. All bytes of the input message are 
converted to ASCII except the first n values as specified 
in the integer input argument. (This was incorporated to 
permit the use of this function with the Motorola GPS 
proprietary message format where the first four bytes are 


gla 


already in ASCII.) 


Input is a pointer of type PANDL and the interger value 
of the number of bytes to skip before conversion. 

Output is a pointer of type PANDL to the new buffer 
with the ASCII converted message. 


ef 
PANDL *Bin_to_ascii( PANDL *bin_msg, int skip ) 
{ 
int i, j, index, tmplen = O; 
ONEBYTE nibble[2], *ptr_tool; 
PANDL *ascii_Msg; 


tmplen = bin _msg->len * 2 - skip; 
ptr tool = bin _msg->ptr; 


/* Create new memory for the converted message. 
sf 

ascil msg = malloc( sizeof( PANDL )); 

ascil msg->ptr = calloc( tmplen, sizeof(ONEBYTE) ); 
ascii msg->len = tmplen; 


/* Copy any initial values to the new buffer that don’t need conversion. 
“af: 
memcpy(ascli_msg->ptr,ptr_tool,skip); 


for (1 = skip; 1 < tmplen; i+ + ) 


/* Split the byte into its upper and lower nibbles. 


eh 
nibble([O] = ( ptr_toolli] > > 4) & OxOF; /* upper */ 
nibble[1] = ( ptr_toolli] & OxOF ); /* lower */ 
for Gj—"O7 = 2. 
{ 
/* Convert the nibbles to ASCII. 
| 
if ( nibble[j] < 10) 
nibble[j] + = Ox30; /* nibble + Ox30 =0x30:0x39 */ 
else 
nibble[j] + = 0x37; /* nibble + 0x37 =0x41:0x46 */ 
ie 
/* Compute the index and insert into the ascii_msg buffer 
yi 


index = 2 * i - skip; 
ascii_msg->ptr[ index |] = nibble[0O]; 
ascil msg->ptr[ index + 1] = nibble[1]; 


}; 


return ascii_ Msg; 


Like 


i 
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Msg_verification() 


Specifically used to verify the integrity of the 
Motorola GPS proprietary message format. 


Input are a pointer of type PANDL and an integer designating 
whether the message was received from the GPS receiver (Rx) 
or is to be transmitted to the GPS receiver (Tx). 

Output is the integer flag indicating whether or not the 
message passed verification. True means it passed. 


a 
{ 
int i, len = -1, array_col; 


int status = TRUE; 


ONEBYTE byteO = msg->ptr[0]; /* Msg prefix @ byte */ 
ONEBYTE bytel = msg->ptr[1]; /* Msg prefix @ byte */ 
ONEBYTE byte2 = msg->ptrl[2]; /* Msg id first byte */ 
ONEBYTE byte3 = msg->ptr[3]; /* Msg id second byte */ 
ONEBYTE byteCS = msg->ptr{ msg->len-3 ]; /* Check Sum byte */ 
ONEBYTE byteCR = msg->ptr{ msg->len-2 ]; /* Carr Retn byte */ 
ONEBYTE byteLF = msg->ptr[ msg->len-1 ]; /* Line feed byte */ 


/* Determine if message verification is for Tx or Rx because the 
number of byte in the message is different for a Tx message 


| 

int Msg_verification( PANDL *msg, int msg_direction } 
| than for an Rx message. 

| 


*/ 
if (msg _ direction ) 
array_col = Q; /* Tx message: to GPS receiver */ 
| else 
array col = 3; /* Rx message: fm GPS receiver */ 


/* Message ID existence verification and expected length. 

wi 
| if (byte2 = 
for (1 


= ‘A’ ) 
O;1 < 30; i+ + ) 


| if ( byte3 = = A_arrayli][array_ coll] ) 


if (msg_direction ) 


len = A_arraylil[1]; /* Msg size to PVT-6 */ 
else 

len = A_arrayliJ[2]; /* Msg size fm PVT-6 */ 
break; 
}; 


Li? 


I 
if (byte2 == ‘B } 
for (1 ="Oe = sidi 
if ( byte3 = = B arraylillarray_ col] ) 


if (msg _ direction ) 


len = B_arrayliJ[1]; /* Msg size to PVT-6 */ 
else 
len = B _arrayliJ[2]; /* Msg size fm PVT-6 */ 
break; 
i 
}; 
iabovie2.——. © 4} 
for) =20) <i ie) 
if ( byte3 = = C_arraylillarray_coll] ) 
{ 
if ( msg _ direction ) 
len = C_array(il[1]; /* Msg size to PVT-6 */ 
else 
len = C_array(il[2]; (P vISGisize Timi say Vea, 
break; | 
iF 
Mi 
if (byte2 != ‘A’ && byte2 != 'B’ && byte2 != 'C’ ) 


{ 
Status = FALSE; 


return Status; 


HE 


/* Message length verification. ("Cj" message is a special case. 
See GPSINIT.H.) 

Ty, 

if (len != msg->len && byte2 != 'C’ && byte3 != ’j’ ) 


status = FALSE; 
return Status; 


he 


/* Message check sum verification. 
ee) 
if (byteCS != Check sum( msg ) ) 


Status = FALSE; 
return status; 


}; 
/* Message begins with @@ prefix. 
mi 
if (byteO != '@’ |! bytel != ‘@’ ) 
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{ 
status = FALSE; 


return Status, 


hi 


/* Message ends with CR-LF. 
oa 
if ( byteCR != OxOD |! byteLF != Ox0OA ) 


status = FALSE; 
return status; 


return Status; 


}; 
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Parse_inbuff() 


This function accepts a pointer of type PANDL to the 
contents extracted from the serial port Rx buffer. As 
this function repeatedly parses each succeeding message 
from this buffer, the buffer’s length decreases until 

the only remaining information is a partial message. 

The pointer is a static pointer and will not lose its 
information upon completion of this function so it will 

be rejoined with the rest of itself when the next grab 
from the Rx buffer occurs. 


Inputs are a pointer of type PANDL and a pointer to a flag 
indicating when the only remaining contents constitute 
a partial message. 

Output is a pointer to a new buffer of type PANDL containing 
the next complete single message in the buffer. Also, 
the inbuff pointer passed into Parse inbuff has been 


altered. 
7, 
PANDL *Parse_inbuff( PANDL *start_of_next_message, int “partial msg _flag ) 
{ 
int i—-O 4 — 70° 
int msg _start_found before = FALSE; 
PANDL *parsed msg; 
ONEBYTE *ptr_tool; 


parsed msg = malloc( sizeof( PANDL ) ); 
ptr tool = start_of next_message-> ptr; 
/* CASE 41: A complete message has been found. The pointer 


tool points to the beginning of a msg (indicated by "@@") 


LZ 


and "@@"” has been found once before. Move the place keeper 
pointer "start of next_msg->ptr” and adjust the length 
"start_of next_msg->len". Copy the “found” message to a 
new buffer and return its location in a pointer of type 

PANDL. 


ei 
for (i = 0; i < start_of_next_message->len; i+ + ) 
{ 
if ( ptr_toolli] = = '@’ && ptr_toolli + 1] == ‘@’) 
{ 
if (msg_start_found before ) 
{ 
parsed msg->len = i-j; 
parsed _msg->ptr = calloc(parsed_msg->len, sizeof(ONEBYTE)); 
memcpy(parsed_msg-> ptr, ptr_ tool, parsed_msg->len); 
Start_of next_message->ptr + = |; 
Start_of_next_message->len -= 1; 
*partial msg _ flag = FALSE; 
return parsed msg; 
msg _ start_found before = TRUE; 
aa 
i 
}; 


/* The start of a second message "@@" was not found. The contents 
are either exactly a complete message or they constitute a 
partial message. Thus, the Parse _inbuff function Is 
given another chance otherwise the "partial-msg_ flag" is 
returned and the partial message contents will be prepended 
to the next buffer brought into Parse _inbuff. 
| 


/* CASE #2: There is exactly one message in the buffer. 
=] 


if ( ptr_tool[0O] == '@' && \ 
ptr_tool[1} == '@’ && \ 
ptr_tooll start_of_ next_message->len - 2} == Ox0D && \ 
ptr_tool{ start_of next_message->len- 1] == Ox0A _) 
{ 
parsed msg->len = i - j; 


parsed msg->ptr = calloc(parsed msg->len, sizeof(ONEBYTE)); 
memcpy(parsed_msg->ptr, ptr_tool, parsed _msg->len); 


start_of_next_message- > ptr 


NULL; 
Start_of next_message->len 


0; 


*partial msg flag = FALSE; 
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/* CASE #3: There is only a partial message left in the buffer. 
a 


else 
{ 
*partial msg flag = TRUE; 
}; 
return parsed msg; 
hi 
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i. 


Master _gps() 


For use with a Motorola GPS receiver configured in the 
master station mode. It builds two binary files on disk: 
One records all output messages and the other records just 
position messages. The function returns a pointer of type 
PANDL that points to a buffer with psuedorange correction 
messages ready for transmission to the slave receiver. 


Input is nothing. 
Output is a pointer of type PANDL. 


PANDL *Master_gps(void) 


{ 
Static PANDL mxsbuff; 
static ONEBYTE mxs[300); /* buffer for partial msg */ 
static ONEBYTE mstart = TRUE; 
ONEBYTE *tail; 
PANDL *inbuff, *tmpbuff, *parse tool, *ben; 
int *partial msg_flag; 


partial msg _ flag = malloc( sizeof( ONEBYTE ) ); 
*partial msg _ flag = TRUE; 


inbuff = malloc( sizeof( PANDL ) ); 
tmpbuff = malloc( sizeof{( PANDL }) ); 
parse tool = malloc( sizeof( PANDL ) ); 
ben = malloc( sizeof( PANDL ) ); 


/* Determine amount of info in Rx buffer and import to tmpbuff. 
ey 

tmpbuff->len = (ONEBYTE) sio_iqueue( MGPS PORT ); 
tmpbuff->ptr = calloc( tmpbuff->len, sizeof(ONEBYTE) ); 

sio read( MGPS_ PORT, tmpbuff-> ptr, tmpbuff->len ); 


/* Initialize the excess buffer that is to contain the partial msg. 
Be 
if ( mstart ) 

{ /* This "if" statement is invoked */ 
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mxsbuff.len = QO; /* only the first time the program */ 
mstart = FALSE: /* is run. Note the scope of mstart. */ 
mxsbuff.ptr = mxs; 


}; 


/* Append tmpbuff messages to the partial message that is in mxsbuff. 
se ) 

inbuff->len = tmpbuff->len + mxsbuff.len; 

inbuff->ptr = calloc( inbuff->len, sizeof(ONEBYTE) ); 

memcpy( inbuff->ptr, mxsbuff.ptr, mxsbuff.len ); 

tail = inbuff->ptr + mxsbuff.len; 

memcpy( tail, tmpbuff-> ptr, tmpbuff->len }); 


/* Free tmpbuff memory and reset partial message buffer length to zero. 
a 

free( tmpbuff->ptr ); 

free( tmpbuff ); 

mxsbuff.ptr 
mxsbuff.len 


mxs; 
0; 


/* Initialize "ben" for use in collecting psuedorange correction messages. 
cy 

ben = inbuff; 

ben->len = 0; 


while ( TRUE ) 
{ 
/* Get a full message from inbuff. 
a) 
parse tool = Parse_inbuff( inbuff, partial _msg_flag ); 
if ( *partial_ msg _ flag ) 
{ 
free( parse tool->ptr ); 
break; 


}; 


/* Verify the received message validity and indicate when 
the message fails verification. 

a 

if (Msg_verification( parse tool, RX )!= TRUE ) 


printf("Message from master receiver fails verification!\n"); 
free( parse _tool->ptr ); 
continue; 


}; 


/* Store only psuedorange corrections in old inbuff buffer for 
return to the calling routine. 
ch 
if (parse _tool->ptr{2] = = ‘'C’ && parse tool->ptr{[3] = = ’e’ ) 
memcpy(ben->ptr + ben->len, parse _tool->ptr, parse tool->len); 
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/* Write all messages to this file. 

so) 

m_all_ file = fopen( “mall.txt", “ab” ); 

fwrite(parse tool->ptr,sizeof(ONEBYTE),parse tool->len,m_all_file); 
fclose( m_all_ file ); 


/* Write only position messages to this file. 
| 
if (parse tool->ptr[2] = = 'B’ && parse tool->ptr[3] = = ‘a’ ) 
{ 
m_pos file = fopen( "mpos.txt", “ab” ); 
fwrite(parse tool->ptr,sizeof(ONEBYTE),parse_ tool- >len,m_pos_file); 
fclose( m_pos file ); 


}; 


free( parse _tool->ptr ); 
}; /* end while */ 


/* Store any partial message in the excess buffer. 
ak 
if ( inbuff->len > O) 
{ 
memcpy(mxsbuff.ptr, inbuff->ptr, inbuff- >len); 
mxsbuff.len = inbuff->len; 


}; 


free( inbuff —); 
free( parse tool ); 
free( partial msg_flag ); 


return ben: 


}; 


ree ee ee ee ee ee ee ee eT ae 
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Slave _gps() 


For use with a Motorola GPS receiver in the slave 
station mode. It converts five data messages into 
scaled and formatted data: 

Position/Status/Data message 

xDOP Table Status message 

Ephemeris Data Output message 

Satellite Range Data Output message 

Position/Status/Data Extension message 
In addition, these and all other output messages are 
returned to the calling routine for transmission to 
the master receiver for post-processing. 


Input is a pointer of type PANDL with any input messages. 
Output is a pointer of type PANDL with all output messages. 
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e 
PANDL *Slave_gps( PANDL *input) 


static PANDL sxsbuff; 
static ONEBYTE sxs[300]; /* buffer for partial msg */ 
static ONEBYTE sstart = TRUE; 
ONEBYTE “*tail; 
PANDL “inbuff, *tmpbuff, *parse tool, *ben; 
int *partial msg _ flag, |; 


/* Deliver any input messages to the GPS receiver. 
a 

sio_ putb( SGPS_ PORT, input->ptr, input->len ); 
free( input->ptr ); 

free( input ); 


/* Create flag pointer and working pointers of type PANDL. 
ro 
partial msg flag = malloc( sizeof(ONEBYTE) }); 


*partial msg flag = TRUE; /* Used in Parse_inbuff */ 
inbuff § = malloc( sizeof( PANDL ) ); 
tmpbuff = malloc( sizeof( PANDL ) ); 


parse_tool = malloc( sizeof( PANDL ) ); 
ben = malloc( sizeof( PANDL ) ); 


/* Determine amount of info in Rx buffer and import to tmpbuff. 
oi 

tmpbuff->len = (ONEBYTE) sio_iqueue( SGPS_PORT ); 
tmpbuff->ptr = calloc( tmpbuff->len, sizeof(ONEBYTE) }); 
sio_read( SGPS_ PORT, tmpbuff->ptr, tmpbuff->len ); 


/* Initialize the excess buffer that is to contain the partial msg. 
a) 


if ( sstart ) 
{ /* This “if” statement is invoked */ 
sxsbuff.len = O: /* only the first time the program */ 
sstart = FALSE; /* is run. Note the scope of sstart.*/ 
sxsbuff.ptr = sxs; 
}: 


/* Append this info in tmpbuff to excess msg from last read. 
a 

inbuff->len = tmpbuff->len + sxsbuff.len; 

inbuff->ptr = calloc( inbuff->len, sizeof(ONEBYTE) ); 
memcpy( inbuff->ptr, sxsbuff.ptr, sxsbuff.len ); 

tail = inbuff->ptr + sxsbuff.len; 

memcpy( tail, tmpbuff-> ptr, tmpbuff->len ); 


/* Free tmpbuff memory and reset partial message buffer length to zero. 
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free( tmpbuff-> ptr ); 
free( tmpbuff }); 
sxsbuff.ptr = sxs; 
sxsbuff.len = O; 


/* Initialize "ben" for use in returning all messages to the calling 


environment. 


i 
ben = inbuff; 
while ( TRUE ) 


/* Get a full message from inbuff. 


a 


parse tool = Parse inbuff{ inbuff, partial_msg_flag ); 
if ( *partial msg_flag ) 


{ 


free( parse tool->ptr ); 


break; 
}; 


/* Verify the received message validity. 


uy 


if (Msg_verification( parse tool, RX )!= TRUE ) 


{ 


free( parse _tool->ptr ); 
continue; 


}; 


/* Convert the binary msg to GPS data structure for use by cpu. 


*/ 


switch (parse_tool-> ptr[3]) 


{ 
Gase a: /* @@Ba Position/Status/Data 
/* Insert conversion here */ 
break; 
Gase) C. /* @@Bc xDOP Table Status 
/* Insert conversion here */ 
break; 
case ‘f’: /* @@Bf Ephemeris Data 
for(i = 5;i < EPH NUM WORDS +5 ; i+ =3 ) 
{ 
hex_ephemeris{iJ{(O] = parse tool->ptr{i + OQ]; 
hex ephemeris[iJ[1] = parse tool->ptr{i + 1]; 
hex_ephemeris{i](2] = parse tool->ptrli + 2]; 
/* Insert conversion here */ 
break; 
Case aq’: /* @@Bg Satellite Range 


/* Insert conversion here */ 


je 


| 


mh 


a 


ue 


break; 


case ‘k’: /* @@Bk Position/Status/Data Ext 
/* Insert conversion here */ 
break; 

default: 
break; 

}; 


free( parse tool->ptr ); 
}; /* end while */ 


/* Store any partial message in the excess buffer. 
ag 
if ( inbuff->len > 0) 
{ 
memcpy(sxsbuff.ptr, inbuff-> ptr, inbuff->len); 
sxsbuff.len = inbuff->len; 
}; 
free( inbuff ); 
free( parse tool ); 
free( partial msg_flag ); 


return ben; 


}; 


#endif /* GPSFUN_H */ 
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G. GPSINIT.H 


/* GPSINIT.H 


This header file contains the Motorola Proprietary 

Binary Message Format characteristcs for GPS 

receivers. The first column contains the second 

letter of the message ID of a message prepared for 

input to the GPS receiver. The fourth column contains 
the second letter of the message ID of a message 
received from the GPS receiver. For all messages sent 

to the GPS receiver except one, there ts a message reply. 

The second column contains the final length of a 
message prepared to be sent to the GPS receiver. The 
third column is the length of a message received from 
the GPS receiver with that particular message ID. 

Note: All message lengths can be expressed with 
one byte except one: The response message Cj ts 294 bytes. 
A unique dummy size equal to 1 byte has been used to 
indicate the length. 

we 


#ifndef —GPSINIT_H 
#define _GPSINIT_H 


#ifndef GPSDEFIN_H 
#include "c:\borlandc\twite\gpsdefin.h" 


#endif 
ONEBYTE A_array[30][4] = {{ a PO we Sa.) 
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1 teow mele. 
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1. Pe Oe oe 
fee lee lee} 
Pee Oe 2 Se ee 
{ 

{ 

{ 

{ 

joeelay on of 42) 
foo 208 Set 
ioe Go. ot), 
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{) 6 2 ie 
{ “Veo Sven 
{> “WIS ea vie 
{ “x Ore eX sey, 
{Ce Wall aaa 
{  “2',11,°11, “zae 
{ “A’, 8, 8 “A .} 
ie B38, 8,8 oe 
eee 9), 9 2 eas 
(DS, Sa bey}, 
ONEBYTE B_array(11](4] ={{ ‘a’, 8, 68, ‘a’ }, 
{-'b, 8, 92a" bie 
1 €,. 8,682, 4c 
{ ‘d’, 8. 23, dase 
{ ‘'e’, 8, 0, "" ae Be->Cb 
1) HieeS)« 280, “tae 
{ (a) 8, 227.0) ie 
{ “Hr, On 24) Eh ee 
Wee eee 
{ ‘f", 80, 0, ’-' },/* Bf->Ce */ 
{ook eB OO tk 3) 


ONEBYTE C array[11][4]= {{ ‘a’, 7, 9, ‘a’ }, 
{ ibe as: 9, ‘h’ }, 
0, 20 ey eedieaye 
aoe 52, ah "k! i, 


{ 

{ 

{ 

{ 

fea oe Ole vil ty 
{ ‘j', 7, 1, ‘j' }, /* Cj resp msg 294 byte */ 
{ 35 Oho 337 Die be Cb ae 

{="-" 0, 52) “e }, yy Ehe = Ce) 

{ ‘~,0, 80, ‘cc })/* Bi-2Cc 7 


#endif /* GPSINIT_H */ 


H. GPSPORTS.H 


jameor or ORTS.A 


Initialize _master_gps() 
Initialize slave _gps() 
Open_ports() 
Close_ports() 


Both initialization functions need to be greatly improved 
to provide an easy to understand and use interface with 
the user. 

ei 


#ifndef  GPSPORTS_H 
#define _GPSPORTS H 


#include <stdio.h> 
#include “c:\borlandc\twite\head-c.h" 
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Initialize master _gps() 


This function presents the current configuration 

settings for the master GPS receiver. The user is 

offered the opportunity to change any of the parameters. 
In addition, the user is requested to indicate which if 

any output files he will want. For those files, the default 
file name is presented and made available for modification 
or change. 


Input is the serial port number through which the master 
GPS communicates. 
Output is nothing. 
mii 


void Initialize master _gps(int mgps) 


{ 

char si[20]={'@’,’@'’,’B’,’a’,1,22,13,10}; 
char s2[20]={'@',’@’,’B’,’h’,1,43,13,10}; 
sio_flush(mgps, 2); 

sio_putb(mgps,s1,8); 

sio_putb(mgps,s2,8); 

printf("Master GPS receiver initialized.\n"); 
return; 


}; 


nO re a ee ee ee ee ee ee a ee ee 


a 


Initialize slave _gps() 


This function presents the current configuration 

settings for the slave GPS receiver. The user is 

offered the opportunity to change any of the parameters. 
In addition, the user is requested to indicate which if 

any output files he will want. For those files, the default 
file name is presented and made available for modification 
or change. 


Input is the serial port number through which the slave 
GPS communicates. 
Output is nothing. 


a 
void Initialize slave gps(int sgps) 


{ 

char s1(20]={’@’,’@’,’B’,’a’,1,22,13, 10}; 
sio_flush(sgps, 2); 

sio_putb(sgps,s1,8); 

printf("Slave GPS receiver initialized.\n"); 
return; 


i 


[OT EEE EE ee ee ee ee a ee 


Open_ports() 


This function initializes the two serial ports 
on the PCD-744 to be used by the master and slave GPS 
receivers. Before initialization, it either loads the 
drivers for the standard and PCD-744 serial ports or 
verifies that they are loaded. It then configures the 
Baud rating, word length, parity, and stop bits as well 
as the modem and flow control. All commands that begin 


with "sio_..." are PCD-744 commands and may be referenced 
in the PCLS-802 PC-ComLIB Serial Communication Programming 
Library. 


The PCD-744 RS-232 8 port intelligent |O card is 
made by ADVANTECH in Sunnyvale, Ca. (408) 245-6678. 

Input is the slave GPS port number. 

Output is nothing. 


ay | 
void Open _ports( int mgps, int sgps ) 
{ 
int BAUD =e oe /* 9600 (p.12 PC-ComLIB) */ 
int IOMODE = (0x03/0x00!0x00); /* 8-N-1 (p. 12) */ 
int MODMODE = Ox00; /* DTR and RTS off (p.26)*/ 
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int HWMODE = Ox0O0O; /* HW and SW flow ctrl off 
(33): / 


n = sio_ioctl{ mgps, BAUD, |OMODE ); 
if (n != QO) 
{ 
printf("Opening comport %d ioctl error.\n",mgps); 
abort(); 
}; 
n = sio_Ictri( mgps, MODMODE ); 
if (n != Q) 


printf("Opening comport %d Ictrl error.\n",mgps); 
abort(); 


n = sio flowctrl( mgps, HWMODE ); 
if (n !}= QO) 


printf("Opening comport %d flowctrl error.\n",mgps); 
abort(); 
i 
n = sio_open( mgps ); 
if (n != Q) 


printf("Opening comport %d error.\n",mgps); 
abort(); 


n = sio_ioctl{ sgps, BAUD, |OMODE ); 
if (n != OQ) 


printf("Opening comport %d ioctl error.\n",sgps); 
abort(); 
he 
n = sio_Ictri( sgos, MODMODE ); 
if (n != O) 
{ 
printf("Opening comport %d Ictrl error.\n",sgps); 
abort(); 
hi 
n = sio_ flowctrl{ sgps, HWMODE ); 
if (n!= QO} 
{ 
printf("Opening comport %d flowctrl error.\n",sgps); 
abort(); 
lr 
n = sio_open( sgps ); 
if (n != Q) 
{ 


ie 


printf("Opening comport %d error.\n",sgps); 
abort(); 
i 
printf("Opened ports %d (Master GPS) and %d (Slave GPS).\n", \ 
mgps, sgps); 
return; 


, 


eee ee ee ee ee ee Se geen eee ee 


Close_ports() 


This function closes and resets all serial ports for 
the PCD-744 Serial |O card. See the Open_ports function 
for more information on the commands used herein. 


Inputs are the master and slave GPS port numbers. 
Output is nothing. 
ak 


void Close ports( int mgps, int sgps ) 
{ 
printf("Closing all ports\n"); 
sio_close( sgps ); 
sio_close( mgps ); 
sio_reset(); 
return; 


}: 


#endif /* GPSPORTS H */ 
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I. 


GPSTRUCT.H 


i GPSTRUCT.H 


ey 


#ifndef 


#define 


#ifndef 


_GPSTRUCT_H 
_GPSTRUCT_H 


_GPSINIT_H 


#include “c:\borlandc\twite\gpsinit.h” 


#endif 


typedef 


struct{ /* Pointer AND Len */ 
ONEBYTE *ptr; 
int len; 


}PANDL; 


/* Struct and global data definitions 


i 


struct T FLOAT EPHEMERIS { 
UNSIGNED _ONEBYTE iode; 


/* Field name abbreviations were extracted from ICD-GPS-200. 


a 


double 


double 
double 
double 
double 
double 


double 
double 
double 
double 
double 
double 
double 
double 
double 
double 


double 


toe; /* Time of ephemeris. si 

moO; /* Mean anomaly at refernce time. a | 

delta _n; /* Mean motion difference from computed value. 

e; /* Eccentricity. a 

sqrt_a; /* Square root of semi-major axis. ss 

omega 0; /* Longitude of ascending node of orbit plane 
at weekly epoch. a 

iO; /* Inclination angle at reference time. a 

WwW; /* Argument of perigee. nh 

omega dot; /* Rate of right ascension. a! 

i_ dot; /* Rate of inclination angle. | 

GUC: /* Amplitude of the cosine harmonic correction 
term to the argument of latitude. ns 

cus; /* Amplitude of the sine harmonic correction 
term to the argument of latitude. | 

CiC; /* Amplitude of the cosine harmonic correction 
term to the orbit radius. oe 

Crs; /* Amplitude of the sine harmonic correction 
term to the orbit radius. i 

CIC; /* Amplitude of the cosine harmonic correction 
term to the angle of inclination. a 

CIS; /* Amplitude of the sine harmonic correction 
term to the angle of inclination. */ 

tgd; /* Estimated group delay differential. we 


f3 5 


ie 


double toc; /* Clock data reference time. | 


double af2; /* Polynomial coefficient (SV clock correction) */ 
double af1; /* Polynomial coefficient (SV clock correction) */ 
double afO; /* Polynomial coefficient (SV clock correction) */ 
double n; /* Corrected mean motion. i 
double e_sart; /* Square root of eccentricity” 2 subtracted 
from 1. = 
double wk1; /* Omega_ dot - WE. eh 
double wkO; /* Omega_O - WE * toe. ss 
double sin_w; /* Sine w. mh 
double cos _w; /* Cosine w. a 
double a; /* Semi-major axis. a) 
double sqrta e fconst; /* Square root of semi-major axis * 
eccentricity * FCONST. ad 


} float_ephemeris ; 


struct T RECEIVER_CHANNELS { 


UNSIGNED _ONEBYTE svid; 
UNSIGNED_ONEBYTE mode; 
UNSIGNED _ONEBYTE strength; 


UNSIGNED_ONEBYTE flags; 
ee 


struct T GEODETIC { 


TWOBYTE degrees ; 
UNSIGNED TWOBYTE minutes ; 


double seconds ; 


ie 


struct T POS CHAN STATUS { 


}; 


UNSIGNED _ONEBYTE month; 
UNSIGNED ONEBYTE day; 
UNSIGNED_TWOBYTE year; 
UNSIGNED _ONEBYTE hours; 
UNSIGNED_ONEBYTE minutes; 


double seconds; 

struct T GEODETIC latitude; 

struct T GEODETIC longitude; 

double datum height; /* meters */ 
double msl height; /* meters */ 
double velocity; /* m/sec */ 
double heading; /* degrees */ 
double current_dop; 


UNSIGNED ONEBYTE dop_ type; 

UNSIGNED ONEBYTE visible_sats; 

UNSIGNED _ONEBYTE sats_tracked; 

struct T RECEIVER CHANNELS channel[NUM_ CHANNELS] ; 
UNSIGNED ONEBYTE rcvr_status; 


#endif /* GPSTRUCT_H */ 


SiG 


J. HEAD-C.H 


/* moxa-c.h ver 1.00 10/11/1991 
* Definitions for MOXA Serial 1/O Controller toctrl 


yh 


/* BAUD rate setting 


#define B50 
#define B75 
#define B110 
#define B134 
#define B150 
#define B300 
#define B600 
#define B1200 
#define B1800 
#define B2400 
#define B4800 
#define B7200 
#define B9600 
#define B19200 
#define B38400 
#define B57600 


/* MODE setting 


#define BIT 5 
#define BIT 6 
#define BIT 7 
#define BIT 8 


#define STOP_1 
#define STOP 2 


#define P_EVEN 
#define P ODD 
#define P_ SPC 
#define P MRK 
#define P NONE 


/* MODEM CONTROL setting 


#define C_DTR 
#define C_RTS 


/* MODEM LINE STATUS 


#define S CTS 
#define S DSR 
#define S RI 
#define S CD 


0x00 
0x01 
0x02 
0x03 
0x04 
0x05 
0x06 
0x07 
0x08 
0x09 
Ox0A 
Ox0B 
0OxO0C 
Ox0D 
OxOE 
OxOF 


-) 
0x00 
0x01 
0x02 
0x03 


0x00 
0x04 


0x18 
0x08 
0x38 
0x28 
0x00 


0x01 
0x02 


Ox01 
0x02 


0x04 


0x08 


/* function declaration 


i) 


/* Word length define 


/* Stop bits define */ 


/* Parity define 


Sok 


ie 
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y 


long sio_iqueue(int); 
long sio_ifree(int); 
long sio_oqueuel(int); 
long sio_ofree(int); 
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APPENDIX F: NPS GEODETIC MARKS 





SPANAGEL GEODETIC MARKS 
WGS-84 


----> 50.85 Deg 


Spanagel Hall runs SW to NE. The marks are on the 8 th 
level near the North-East cornor of the building. 


Marks 1 and 2 are standard survey disks, marked SPA-1 and 
SPA-2. SPA-3 is a square of thin aluminum. All marks 
are attached to the roof under a rail. 


Poles are mounted on the r&il above the marks. The 
height to the marks are 11.83 m ( WGS-84 elipsoid ). 


Mark Pole 
SPA-1 36 35.70375 121 52.47777 11.83 14.08 
SPA-2 36 35.70180 121 52.48074 11.83 14.10 


SPA-3 36 35.70277 121 52.47926 11.83 14.09 


BESS, 


Table IX WGS 84 SPANAGAL HALL EIGHTH FLOOR GEODETIC MARKS. 







SO — 













Loca~ Item degrees degrees degrees meters 
minutes minutes (above 
seconds geoid) 


j Lat | 36.59506253 136° 35.70375'| 36° 35° 
| °N Oa ay ed 
j Lon {121.87462956 A WA ays, 
°'W 52.47777° 28.6664”’ 

Lat | 36.59503006 |36° 35.70180"} 36° 35’ 
°N 42.1082”’ 

Lon |121.87467903 a) 2 eS 8 
”'W 52.48074’ 28.8445”’ 


| 11.83 
spa-3| Lat | 36.59504622 [36° 36.70277 36° 35’ 
°N 42.1664" 

2323 | Lon |121.87465439] 121° 121° 52° 

OW 52.47926" | 28.7558” 
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